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ABSTRACT 


Intensity  interferometry,  as  developed  by  Banbury  Brown  and  Twiss 
for  stellar  observation,  has  shown  relative  insensitivity  to  atmospheric 
scintillation.  However,  with  classical  sources,  the  limitations  placed 
on  this  technique  by  quantum  noise  r.r.d  detector  efficiency  arc  severe. 

This  situation  is  vastly  improved  when  laser  illumination  is  employed. 

A  modification  of  the  Hanbury  Brown  -  Twiss  experiment  is  described  whereby 
intensity  correlat iv*i  is  performed  using  coherent  light  and  the  statis¬ 
tical  averages  arc  taken  in  the  space  domain.  Generalizing  a  form  of 
the  mutual  coherence  function,  the  far-zone  behavior  of  the  mutual  inten¬ 
sity  function  for  an  intermediate  time  average  is  derived.  This  result 
is  used  to  reconstruct  the  irradiance  distribution  of  a  spatially  rough 
source.  The  far-field  intensity  distribution  is  recorded  spatially  for 
one  time-resoluticn  unit  of  the  detector.  Tne  resulting  spatial  signal 
is  autocorreiated  and  related  to  the  intensity  distribution  over  the 
source.  Thus,  without  averaging  in  the  time  domain,  a  spatial  Fourier - 
transform  relation  is  derived  oetween  the  far-field  intensity  correlation 
and  the  source  irradiance,  similar  to  ♦'he  results  of  Hanbury  Brown  and 
Twiss. 

A  major  limitation  of  intensity  interferometry  in  both  the  space 
and  time  domains  has  been  that  only  sources  of  even  symmetry  could  be 
uniquely  inferred  from  far-field  intensity  correlations,  5ir.ee  only  the 
modulus  was  derived.  However,  this  restriction  can  be  reisved  by  pro¬ 
cessing,  in  addition,  the  field  in  a  simple  way  before  detection  to 
a  new  symmetrized  function.  This  intensity  record  corresponds  tc  the 
radiation  pattern  derived  from  the  pure  even  part  of  the  intensity  pro¬ 
file.  The  record  is  autocorrc lated ,  as  befora,  yielding  a  spatial  power 
spectrum.  But  by  the  central -ordinate  and  herraitisn  properties  of  the 
spatial  Fourier  transform,  the  real  part  of  the  transformation  can  be 
specified  exactly  from  the  symmetrized  record.  Using  the  modulus  of  the 
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total  transform  together  with  the  real  part  of  the  transforms,  one  can 
infer  the  phase  to  within  a  sign.  Thus  a  pair  of  images  is  derived,  one 
erect,  the  other  inverted.  This  ambiguity  can  be  resolved,  however,  by 
translating  the  source  in  a  known  direction  off  axis. 

The  results  derived  for  the  spatial  defection  of  intensity  fluctu¬ 
ations  are  shown  to  be  similar  to  previous  results  dealing  with  the  phe¬ 
nomenon  of  laser  speckle  in  terms  of  the  spatial  transform  relations. 
However,  the  criterion  of  dtcector  time  resolution,  central  to  the  resulcs 
of  intensity  interferometry,  is  absent  for  laser  speckle  effects.  Thus, 
it  can  be  inferred  that  the  utility  of  rntensity  interferometry  is  based 
on  the  detection  of  a  time-variant  speckle  pattern. 

Finally,  experimental  evidence  is  offered  to  support  the  spatial 
transform  relations  predicted  here.  In  addition,  a  time-domain  intensity- 
correlation  experiment  it  performed,  the  purpose  of  which  is  to  infer  the 
existence  of  a  time-varying  speckle  pattern. 
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PREFACE 


The  history  of  intensity  interferometry  is  rooted  in  the  work  of 
Hanbury  Brown  and  Twiss.  Their  earliest  investigations*  dealt  with  the 
problem  of  resolving  stellar  radio  sources  by  a  technique  involving  the 
correlation  of  the  squared  outputs  of  two  receivers.  The  advantages  are 
the  reduction  of  certain  kinds  of  experimental  constraints  as  well  as  2 
the  comparative  insensitivity  of  the  method  to  atmospheric  scintillat  :on. 

A  preliminary  conclusion  reached  at  that  time  was  that  this  technique  of 
intensity  correlation  would  not  be  applicable  at  optical  frequencies  ^ 
because  of  limitations  imposed  by  photon  noise.  However,  in  later  work, 
Hanbury  Brown  and  Twiss  showed  that  meaningful  intensity  correlations 
could  be  made  at  optical  frequencies  even  with  highly  degenerate  sources. 
The  limitations  placed  on  this  approach  by  quantum  noise  and  detector 
efficiency  have  been  severe,  calling  for  highly  refined  experimental 
technique. 

For  laser  illumination,  the  situation  is  very  different.  The  signal- 
to-noise  ratio  can  be  typically  increased  by  six  orders  of  magnitude.^ 
However,  the  statistics  of  the  source  must  be  considered  in  the  measure¬ 
ment.  The  key  to  relating  intensity  correlations  to  some  property 
involving  field  correlations  lies  in  the  assumption  of  gaussian  statis¬ 
tics,-’  for  which  all  higher  moments  are  determined  from  che  first  and 
second.  Single-mode  lasers,  though,  are  distinctly  non-gaussian  in  their 
temporal  statistics  and,  therefore,  cannot  be  described  by  theory  framed 
for  thermal  sources.  But  with  the  addition  of  axial  modes.,  it  has  been 
asserted6  that  the  field  amplitude  becomes  nearly  gaussian  distributed. 

The  principal  formula  used  by  Hanbury  Brown  and  Twiss  [Ref.  3b, 

Eq.  (2.1)]  to  infer  the  diameter  of  a  distant  source  shows  that  the  time- 
averaged  correlation  of  intensities  at  two  points  is  equal  to  the  product 
of  a  function  involving  the  temporal  characteristics  of  the  source  with 
the  square  of  the  spatial -Fourier  transform  of  the  source  intensity  dis¬ 
tribution.  Consideration  of  the  source  temporal  statistics  is  necessary 
if  the  intensity-product  output  of  the  detectors  is  averaged  in  the  time 
domain  (as  it  nearly  always  is)  to  overcome  the  limitations  imposed  by 
photon  and  detector  noise  and  possibly  reduce  scintillation  effects  pro¬ 
duced  by  transmission  through  the  atmosphere.  If  the  intent  of  an  inten¬ 
sity  correlation  experiment  is  to  gain  information  concerning  the  source 
intensity  distribution,  then  the  temporal  statistics  may  be  of  little 
interest  in  themselves. 

Intensity  interferometry  can  be  understood  as  a  two-point  correlation 
of  intensities  following  the  squaring  of  the  electric  field  at  the 
detector.  If  the  source  is  quasi -monochromatic,  each  differential  ele¬ 
ment  on  the  object  emits  a  number  of  temporal  modes  that  interfere  with 
each  other  at  the  detector.  If  the  detector  has  sufficient  bandwidth  to 


*  References  cere  listed  on  page  66. 
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detect  the'--'  beat  frequencies,  the  amplitude  and  phase  of  the  incoming 
intensity  beats  are  utilized.  What  is  sufficient  depends  on  the  band¬ 
width  of  the  source.  For  a  thermal  source,  most  beat  frequencies  are 
too  high  to  be  resolved  even  with  the  megacycle  response  of  the  Hanbury 
Brown  -  Twiss  apparatus.  For  a  single-m^s  laser,  all  beat  frequencies 
could  be  1  ss  than  100  Hz.  It  is  com\ionly  argued  that,  in  these  circum¬ 
stances,  t  i  random  fluctuations  of  ti  e  temporal  statistics  from  dif¬ 
ferent  poi  .ts  of  the  source  cause  the  beat  frequencies  from  each  source 
point  to  add  incoherently  at  the  detector.  The  time-averaged  intensity 
correlation  is  then  proportional  to  the  squared  spatial  Fourier  trans¬ 
form  of  the  source  intensity.  This  approach  gives  essentially  a  squared 
version  of  the  van  Cittert-Zemike  theorem.'7 


We  wish  to  suggest  that  the  requirement  of  surface  roughness  at  the 
source  (to  assure  spatial  incoherence)  is  sufficient  to  guarantee  the 
incoherent  addition  of  beat  frequencies  at  the  intensity  detectors. 

Thus,  if  temporal  noise  (photon  noise,  time -dependent  detector  noise) 
is  largely  absent  in  a  local  spatial  sense,  as  might  be  the  case  with  a 
multi-axial-mode  laser  with  photographic  detection,  then  the  intensity 
information  might  be  gathered  during  one  resolution  time  of  the  detector 
over  a  plane  section  normal  to  the  direction  of  light  propagation.  Any 
noise  arising  in  the  process  would  be  spatial,  and  might  be  averaged  out 
by  taking  a  sufficiently  large  area  of  spatial  correlation.  The  reduc¬ 
tion  of  atmospheric  spatial  noise  would  be  similar  to  a  process  known  as 
aperture  averaging.®  Film-grain  noise  would  be  extremely  well  averaged 
by  the  relatively  large  area  of  averaging. 


The  relative  insensitivity  of  intensity  interferometry  to  turbulence 
can  be  shown2  if  the  propagation  medium  is  assumed  to  be  dispersionless 
over  the  range  of  carrier  frequencies.  Because  each  temporal  frequency 
sees  the  same  refractive  index,  the  differential  (beat)  frequencies 
remain  unchanged.  However,  the  spatial-Fourier-transform  relation  between 
the  source  and  the  far-field  scales  as  the  average  frequency,  not  the 
beat  frequency,  and  thus  the  resolution  afforded  by  optical  frequencies 
is  maintained. 


The  idea  of  examining  spatial  beat  frequencies  of  second-order 
correlation  is,  of  course,  not  new.  Many  classical  field-correlation 
interferometers,  as  well  as  holographic  experiments,  are  built  on  this 
principle,  involving  a  spatial  or  time  lag  between  interfering  beams  of 
the  same  source.  More  difficult  is  the  spatial  recording  of  beats  from 
two  independent  sources,  as  demonstrated  by  Magyar  and  Mandel . 


In  this  thesis,  we  examine  the  subject  of  intensity  interferometry 
in  the  spatial  domain.  In  Section  1,  an  introduction  to  this  field  is 
given,  leading  up  to  the  results  of  Hanbury  Brown  and  Twiss  and  the 
extension  of  their  work  by  Gamo.  In  Section  2.  using  as  a  basis  a 
mathematical  framework  due  to  Marchand  and  Wolf,^-  we  generalize  the 
concept  of  the  mu' •  1  ‘-'■'herence  function  to  time  averaging  of  arbitrary 
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length.  The  fourth-order  field  correlation  function  in  the  far  zone  is 
then  formed,  revealing  a  Fourier  transform  relationship  with  the  source 
irradiance. 

In  Section  3,  a  technique  of  symmetrizing  the  electric  field  before 
detection  is  developed  as  a  way  of  inferring  the  phase  of  the  source- 
intensity  Fourier  transform.  With  this  additional  information,  the 
transform  can  be  inverted  to  form  the  intensity  distribution  for  sources 
of  arbitrary  symmetry.  In  the  next  section,  comparisons  are  made  between 
spatial  intensity  interferometry  and  the  phenomenon  of  laser  speckle. 
Finally  in  Section  5,  some  experiments  involving  spatial  intensity  cor- 
relati -n  are  described.  Included  is  an  investigation  into  the  relation¬ 
ship  between  the  average-value  and  time-varying  speckle  patterns. 


1.  AN  INTRODUCTION  TO  INTENSITY  INTERFEROMETRY 


1.1  Introduction 


Central  to  the  framework  of  intensity  interferometry  are  two  kinds 
of  relationships,  each  existing  in  distinct  domains.  The  first,  and  the 
principal,  of  these  relationships  lies  in  the  spatial  domain,  since  we 
are  here  primarily  interested  in  optical  imaging  by  means  of  intensity 
correlations.  Specifically,  by  some  manipulation  of  the  far-field  inten¬ 
sity  distribution  of  a  source,  we  wish  to  infer  the  intensity  distri¬ 
bution  over  the  source  itself.  However,  the  efficacy  of  this  intent 
rests  almost  totally  on  the  constraints  of  a  second  area  of  relationships, 
the  temporal  domain.  It  is  here  that  the  bandwidth,  power,  and  time 
stability  of  the  carrier  waves  are  determined  and  even  the  practicality 
of  this  special  kind  of  imaging. 

In  this  chapter,  we  intend  to  give  a  systematic  introduction  to 
intensity  interferometry.  In  essence  this  technique  is  straightforward. 
The  fundamental  principles  are  rc.-ealed  in  the  well-known  mathematical 
techniques  of  correlation  and  Fourier  transformation.  First,  we  will 
establish  the  basic  spatial  relationships  between  source  and  far  field 
for  both  the  first-  and  second-order  cases.  Next,  we  will  show  the 
relationship  between  second-  and  fourth-order  field  statistics  for  a 
gaussian  time  variable.  Finally,  we  will  describe  the  work  of  Hanbury 
Brown  and  Twiss  and  the  extension  of  their  work  by  Gamo. 

1.2  The  Huygens -Fresnel  Principle 

We  wish  to  examine  here  the  relationship  between  the  electric  field 
on  a  section  of  plane  normal  to  the  direction  of  energy  transport  with 
the  field  in  a  plane  section  at  some  later  time.  Referring  to  Figure  1.1, 
we  call  V(.£,n)  the  electr  ic  field  bounued  by  the  aperture  E  in  the  £-n 


Figure  1.1  The  Coordinate  Axes  fc  the  Object  (c-n) 
and  Detection  (x  y)  F  -nes 
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plane.  Upon  realizing  that  each  point  in  the  apeiture  radiates  a  spheri¬ 
cal  wave  to  the  right,  we  can  directly  write  the  Huvgens-Fresnel  principle 
(Ref.  7,  Ch.  8)  for  paraxial  waves  for  which 


V(x,x;t)  -Jj 


V(C.n;t) 


exp  li 


[ikKU, n;x,y, 
iAR(£,n;x,y) 


df.dn, 


(1.1) 


where  k  =  2n/A,  and  a  equals  the  wavelength  of  the  radiation.  Equation 
(1.1)  states  that  the  electric  field  V(x,/;t)  is  formed  of  a  superposi¬ 
tion  of  waves  emanating  from  each  point  within  the  aperture,  properly 
phase-shifted  according  to  the  exponential  term  and  diluted  by  the  1/R 
expression.  Now 


R(e,n;x,y)  =  [Rq4  ♦  (x  -  fj  +  (y-n)2]' 


(1.2a) 


(1.2b) 


when  approximated  by  the  first  two  terms  of  the  binomial  (the  Fresnel 
approximation)  expansion.  Using  Equation  (1.2b)  in  Equation  fl.l)  and 
making  the  far-field  (Fraunhofer)  approximation  [Rn  »  k(£2  +  r,  )  /2] , 

we  can  write 


V  (x,y;t) 


exp[ ik (Rq  (x2  ♦  y 2)/2Rq) ] 


no 

jj  V(5 ,  n;t)exp£-i  (x?+ yn)Jd^dn. 


(1.3) 


Apart  from  the  coefficients  before  the  integral,  Equation  (1.3)  shows 
that  the  electric  fields  in  an  aperture  and  the  far  field  are  related 
by  a  spatial  Fourier  transform  operation.  This  property  is  basic  to 
the  ensuing  work  involving  intensity  correlations  with  quasimonochromatic 
light. 

1.3  The  van  Cittert-Zernike  Theorem 

Although  Equation  (1.3)  holds  for  radiation  of  general  frequencies, 
the  relationship  it  expresses  is  somewhat  academic  here  from  the  stand¬ 
point  that  it  describes  electric-field  quantities  that  are  unmeasurable 
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at  optical  wavelengths.  It  is  this  constraint  that  has  prompted  many 
investigators  to  couch  optical  theory  in  the  form  of  correlations  of 
field  quantities.  The  most  well  known  of  these  correlationr,  is  ce  u  i 
the  mutual  coherence  function,  r(x  ,x  >T).  where 


r(x  ,x  ,t) 
—1  —2 


lim 

T-k» 


1 

2T 


A- 

-T 


t*t)V 


(x2,t)dt. 


(1.4a) 


=  Vfx^.t  ♦  x)V*(x^,t),  (1.4b) 

and  *i  is  the  time  delay  between  the  instantaneous  product  of  the  elec¬ 
tric  fields  at  the  points  Xj  and  x^.  A  special  case  of  the  mutual 
coherence  function  results  when  the  time  delay  is  set  to  zero  or 

r(i.1»£2,1^i=0  E  ^ j  ’  — 2 ^  ’  (1-5) 

where  JQc^.x^)  is  called  the  mutual  intensity  function. 

We  wish  to  calculate  the  mutual  intensity  in  the  far  field  of  a 
spatially  incoherent  source.  Using  Equations  (1.5),  (1.4b),  and  (1.3), 
we  write 


V(x1,yl;t)V*(x2,y2;t)  =  J(x} ,yj ;x2 ,y2) 


(1.6a) 


exp 


fik  2  2  2  ?  I 

hr-  x 2 *  y  1  - 


a> 

//// 


V C £ j  > pj »t)V  (C2,n2;t) 


exp 


Li  jl  „ 

L  Ro 


(1.6b) 


^1  ~  X2^2  *  -vini -y2n2) 


)  d^d52dnjdn2. 


where  the  processes  of  temporal  averaging  and  spatial  integration  have 

been  interchanged  and  the  radiation  assumed  quasimonochromatic 

^Vmin  >=>  vmax  "  Vmin^  so  t*iat  t*ie  wave^en8t^1  dependence  can  be  approxi 

mated  by  the  mean  wavelength.  Since  the  source  is  spatially  incoherent 
the  mutual  intensity  takes  the  form1^ 

VU1,n1;t)V*(52,u2;t)  =  J(5j ,  ;S2>  n2) ,  (1.7a) 


IU,n)6U.  -  i  '«(n.  -  nj. 


(1.7b) 


Physically,  Equation  (1.7b)  implies  that  the  time  fluctuations  of  the 
electric  fields  at  two  non-identical  points  in  the  source  plane  are  com¬ 
pletely  uncorrelated;  equivalently,  the  total  power  measured  at  a  point 
in  the  far  field  is  simply  the  sum  of  the  squared  electric  fields  from 
each  differential  element  of  the  source,  taken  with  the  proper  phase 
delay  2no  attenuation. 

Using  £c'»ai:o;»  (1.7b)  in  Equation  (1.6),  the  mutual  intensity  col¬ 
lapses  to  a  single  ar^a  integral  giving 


(1.8) 

*  JJ  I(£,r)exp  j -i  ~  f(Xj  -  x2)£  ♦  (y1  -  y2)n]jd£dn. 

(  O  ) 


JF(x1,y1;x2,y2)  = 


f.  k  ,  2  2  2  2  " 

“p[>  iro  tx,  -  x2  *  n  •  ^ 


(X\)2 


Although  usually  presented  in  a  normalized  form  (Ref,  7,  p,  510), 

Equation  (1.8)  is  a  statement  of  the  van  Cittert- Zernike  theorem  which, 
aside  from  the  accompanying  coefficient,  shows  that  the  mutual  intensity 
in  the  far  fielc  f  a  spatially  incoherent  source  is  given  by  the  Fourier 
transform  of  the  intensity  across  that  source.  This  theorem  is  basic  to 
most  optical  imaging,  since  lenses  effect  the  inverse  transformation  of 
Equation  (1.8)  to  give  a  scaled  distribution  >jf  the  object  irrsdiance. 
fiy  the  ..caling  property  from  one  domain  to  another,  intrinsic  tc  the 
Fcurer  transform  operation,  large  detail  in  the  source  is  given  by  small 
separations  in  the  far  field  and  vice  versa.  Thus  imaging  with  a  finite 
aperture  implies  a  finite  limit  to  the  high-frequency  detail  resolvable 
on  the  source.  This  justifies  the  well-known  description  of  lenses  as 
low-pass  filters.  Also,  since  there  are  many  more  pairs  of  points  within 
an  aperture  corresponding  to  small  separations  than  to  large,  there  is  a 
built-in  redundancy  weighted  in  favor  of  low-frequency  resolution. 

Equation  (1.8),  although  a  statement  of  second-order  correlation, 
is  basic  also  to  fourth-order  correlation,  as  we  shali  see  later. 

1.4  A  Theorem  Concerning  the  Fourth-Order  Gaussian  Rand--  Process 

We  have  just  seen  in  Equation  (1.7)  the  way  in  which  a  time  average 
can  be  usea  to  impose  a  condition  cn  the  correlation  of  an  electric-field 
pair.  No  particular  assumption  was  made  about  the  statistics  of  the 
field  variables.  However,  it  is  well  known  that  classical  thermal  sources 
exhibit  statistical  fluctuations  that  are  gsussian  in  nature.  Hodara6 
has  asserted  that  lasers  with  but  a  few  axial  modes  are,  to  a  good 
approximation,  gaussian  as  will .  However,  Troup13  has  argued  that  gaus- 
sian  statistics  are  achieved  only  in  the  limit  of  a  large  nW-er  of 
axial  modes. 


16 


A  well-known  property  of  gaussian  statistics  is  that  all  higher  order 
moments  are  representable  in  terms  of  the  first  and  second.  A  relation 
given  by  Middleton^  reflecting  this  property  gives  the  expectation  value 
of  the  2mth-order  correlaf’on  for  the  zero-mean  random  variable  zj  = Xj  -  xj, 
where 


V  z  [,v =*]■ 

all 

(1.9a) 

pairs 

and  E(z,z,. . ,z.  , )  =  0. 

i  2  2m+l 

(1.9b) 

h'e  now  use  Equation  (i.9a)  to  examine  the  relation  between  the  fourth- 
order  and  lower-order  moments  of  a  comp lex -valued  electric  field  variable 
that  is  gaussian  in  the  domain  in  which  the  averaging  takes  place.  We 
therefore  write 

E(I1I2)  =  E(V1VIV2V2) 

(1.10a) 

=  V,VT  V2v*  .  V,V2  V'V'  .  v‘v; 

(1.10b) 

*  |V,V2|2  .  Iv,v* !2. 

(1.10c) 

where  the  complex  conjugation  has  been  used  to  assure  a  real  expectation 
value.  The  first  term  of  Equation  (1.10c)  is  the  product  of  the  mean 
intensities  at  points  one  and  two,  while  th»  last  is  the  mutual  coherence 
function  defined  by  Equation  (1.4a)  if  the  overbars  are  taken  as  an 
average  in  the  time  domain  and  the  ergodic  hypothesis  is  assumed.  We 
now  examine  the  second  term  in  this  expression. 

The  electric  field  can  be  represented  by  the  complex 
(where  a  constant  term  has  been  suppressed)  as 

analytic  signal 

rm 

V(t)  =  /  v(u)  exp(-iut)da), 

*'n 

(1.11) 

0 


where  v(w)  exists  only  for  positive  frequencies.  Using  Equation  (1.11), 
the  second  term  of  Equation  (1.10c)  can  be  expressed 
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Jf  J  jj  v1('“Pv2(w2)cxp[-i(a)l  *  u>2)t]dti»i 


-T  0 


05 

£  !//'  ■• 


4  du  dt  ,  (1.12a) 


<-.j)v2(u2)sincf  (uj  ♦  u2)T]dUldu2|  , 


(1.12b) 


where  in  going  from  Equation  (1.12a)  to  (1.12b)  the  order  of  t.jne  and 
frequency  integrations  has  been  reversed,  and  sine  x  =  sin  x/x.  As  the 
limit  of  T  is  taken,  the  sine  function  approaches  a  delta  function  such 
that  the  integral  has  value  only  for  ±Wj  =  +u>2.  Since  v(u)  can  exist 
only  for  positive  frequencies,  the  integral  is  zero. 


corre 


Thus  the  fourth-order  correlation  of  electric  fields  (second -order 
lation  intensities)  given  in  Equation  (1.10a)  becomes 


1 1 1 2  ~  ll  *2  *  *ri?  I  * 

where  the  iefin:ng  relation  of  Equation  (1.4a)  has  been  used. 


(1.13) 


Equation  (1.13),  therefore,  describes  the  relationship  between  the 
intensity  and  the  field  correlations  for  a  process  that  is  gaussian  in 
the  time  domain.  We  note  that,  in  general,  the  field  correlation  is  a 
complex  quantity,  so  that  only  the  relationship  between  the  intensity 
correlation  and  the  modulus  of  the  field  correlation  is  implied.  It 
reveals  the  underlying  principle  by  which  intensity  correlations  in  the 
far  field  may  be  used  to  infer  the  accompanying  field  correlations  and 
hence,  through  Equation  (1.8),  to  gain  knowledge  of  the  intensity  dis¬ 
tribution  at  the  source. 

1.5  Intensity  Interferometry  in  the  Temporal  Domain 

As  indicated  initially  in  the  Preface,  Hanbury  Brown  and  Twiss 
were  the  first  to  realize  the  potential  of  using  intensity  beats  to 
derive  information  concerning  the  irradiance  of  a  distant  source.  First 
working  at  radio  frequencies  and  then  later  at  optical  wavelengths. ^ 
Hanbury  Brown  and  Twiss  formulated  the  imaging  problem  and  used  the 
results  to  infer  the  diameter  of  the  star  Sirius.  We  wish  to  present  a 
brief  outline  of  their  work  which  relates  to  the  problem  at  hand. 


Following  Reference  3b,  Hanbury  Brown  and  Twiss  express  the  elec¬ 
tric  field  as  a  superposition  of  quaslmonochromatic  waves  emitted  from 
a  diffeiential  area  in  the  form 


V  (x)  « 


E 


(1.14) 


2 

where  hr(^)d^  is  proportional  to  the  population  of  the  rth  mode,  c  is 
the  vacuum  velocity  of  light,  T  is  the  time  of  observation,  and  x^.(£)  is 
a  random  phase  that  is  uniformly  distributed  between  0  and  radians. 

This  phase  is  both  time  and  space  dependent.  The  condition  of  spatiai 
incoherence,  expressed  earlier  in  the  form  of  Equation  (1.7),  was  assured 
by  Hanbury  Brown  and  Twiss  by  requiring  that 

xrCOxs(i’)  *  Ci-is)  j 

i 

Using  Equation  (1.14),  the  instantaneous  intensity  is  formed,  where  1 


OD 

(i.t)  *  j j  *iti!  g  | J  / 


2e(alralsnlrrils)iS 


f 1 . 16) 


*  cos 


-  ff{rR(£,x)  -  sR(£,x)}  -  |xr(£)  -  Xs(£' ) |]  , 


nlr  is  proportional  to  the  number  of  quanta  per  mode,  alr  is  the  detector 
quantum  efficiency,  and  e  is  the  electronic  charge;  also,  following 
Hanbury  Brown  and  Twiss,  two-dimens i^nal  spatial  integration  is  indicated 
by  single  integral  signs.  Ii  arriving  at  Equation  (1.16),  the  form  of 
Equation  (1.14)  was  square!,  giving  a  sura  of  two  cosine  terms,  one 
involving  a  sum  frequency,  the  other  a  difference  frequency.  The  former 
term  has  been  dropped,  since  it  corresponds  to  a  very  high  frequency 
term  (approximately  twice  the  mean  optical  frequency).  The  second  term 
can  be  visualized  in  the  form  of  a  square  sxr  symmetrical  matrix.  The 
diagonal  (r=s)  corresponds  to  a  dc  term  that  has  also  been  dropped,  since 
che  ensuing  electronic  apparatus  has  zero  dc  response.  Twice  the  sum  is 
taken  of  the  off-diagonal  terms  in  the  upper  half  of  the  matrix. 

Next,  to  account  for  the  temporal  resolution  of  the  apparatus,  the 

complex  freauency  response  is  represented  by  so  that  the  filtered 

current  J(x,t)  is  written  ' 
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J(x.t) 


Off 

■  h dyi'  e  ld- 

J  z  r>s  s= l J 


^e^°lralsnIrnls^ 


Invoking  the  Fresnel  approximation  [Equation  (1.2b)]  for  the  ray  paths 
R(E,x),  Equation  (1.17)  is  used  to  form  tne  t'me-average  cross-correlation 
function 


where 


(1.18) 


=  b.  QLi  -  *•>).  f '  =  *  "  * 


(1.19) 


and  k  is  a  special  kind  of  vector  wave  number  defined  by  Hanbury  Brown 
and  Twiss  having  dimensions  of  time/length.  Under  the  assumption  of 
ergodicity,  the  average  defined  by  Equation  (1.18)  is  taken  as  the  equiva 
lent  of  the  ensemble  average. 

When  Equation  (1.17)  is  used  to  form  the  product  of  Equation  (1.18), 
a  pair  of  cosine  terms  result,  similar  to  the  step  from  Equation  (1.14) 
to  Equation  (1.16).  The  first  cosine  argument  involves  a  sum  of  phase 
terms,  the  second  a  difference.  Reflecting  on  the  condition  expressed 
by  Equation  (1.15),  it  becomes  apparent  that  only  cosine  terms  inde¬ 
pendent  of  the  random  phases  can  contribute  to  the  integral  expressed 
by  Equation  (1,18).  The  cosine  composed  of  sum  terms  therefore  drops 
out,  leaving  the  latter  term  with  the  additional  constraint  of  Equation 
(1.15)  forcing  an  equivalence  among  pairs  of  frequency  variables,  simi¬ 
lar  to  the  earlier  result  expressed  by  Equation  (1.7). 

The  result  of  the  operation  defined  by  Equation  (1.18)  is 
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to 

cjp)  »  dWi'  dxdx '  J"  dvai(v}«2(v) 


x  [Hj  <r.v,On2(v1,On1  (v,£')n2(v, £_*)]' 


(1.20) 


i  2v 

cos  .  “tt 
( cR. 


-CO  • 

(x-x')  (£-£))•  j  ^{F1(f)Fj(f)  +  F*(f)F2(f)| 


a  relationship  that  can  more  succinctly  be  expressed  in  the  form 

ccTT  =  |r(v,f‘)|2  cW  ,  (1.21) 


irtv.f')!2  .-~- 


ff*& i’  I(OKc’)  COS  j~~i'  •  (£-£’)| 

cx> 

ffdk6L'  1  (£)!(£) 


»  (1.22) 


and  C(0)  is  the  correlation  at  zero  spacing  of  the  two  intensity  detec¬ 
tors.  Equation  (1.21)  is  the  principal  result,  revealing  a  particularly 
important  relationship.  It  is  shewn  that  the  process  of  intensity  corre¬ 
lation  in  the  far  field  with  an  infinite  time  average  can  be  formulated 
as  a  product  of  two  distinct  terms.  The  first  involves  the  square  of  the 
Fourier  transform  of  the  irradiance  distribution  across  the  source.  The 
second  term  assumes  the  role  of  proportionality  constant  and  relates  to 
the  power  of  the  beat-frequency  carrier. 

_  We  note  that  the  source  spatial  information  and  coefficient  term, 

C(0),  are  completely  decoupled  for  the  above  condition  of  temporal 
averaging.  This  important  result  leads  us  to  speculate  under  what  con¬ 
ditions  an  illuminated  rough  surface  might  enjoy  a  similar  description 
of  decoupled  space-time  statistics  for  the  case  in  which  short  time 
averaging  takes  place.  In  Section  2.  this  subject  will  be  explored. 

Finally,  we  see  that  Equation  (1.22)  gives  the  squared  modulus  cf 
the  spatial  Fourier  transform  of  the  source  distribution.  This  means 
that,  unfortunately,  except  for  certain  situations  for  which  the  phase 
of  the  transform  is  known  a  priori.  Equation  (1.22)  cannot  be  inverted 
to  find  I(£).  A  method  to  circumvent  this  constraint  is  the  subject  of 
the  next  section. 


1 . 6  The  Triple  Correlator  of  Gamo 

As  have  found  in  the  two  previous  sections,  fourth-order  field 
correlation  yields  information  only  about  the  absolute  square  of  the 
second-order  correlation.  Equivalently,  only  a  power-spectral  measure¬ 
ment  is  made  relating  to  a  function  that  is,  in  general,  complex.  Thus 
it  -a  impossible  to  invert  the  power  spectrum  to  derive  the  far-field 
mutual-intensity  function  and  thence  the  irradiance  distribution  over  a 
spatially  incoherent  surface. 


In  order  to  gain  information  about  the  phase  of  the  mutual  coherence 
function  by  the  method  of  intensity  correlations,  Gamo^  proposed  a 
sixth-order  correlation  technique.  We  discuss  briefly  the  principle. 
Using  Equations  (1.9a)  and  (1.4a),  the  third-order  correlation  of  inten¬ 
sities  becomes 

Edj^Ij)  =  ^  tv  j v  * v  2 V  2V  3V  3 )  (1.23a) 


■  T3.T(|r12i*.  |r23l2  *  |r31|2) 

*  2lr15l  i  r23  !  Ir3!*  COS(*l2**23**31>'  (1  •  23b) 


whore  I,  =  1=1=  I,  and  the  terms  in  <f., , 
12  3  1} 

between  points  i  and  j.  Using  the  identity 


indicate  the  phase  difference 


Al j AljAlj 


■  VA  •  ‘3  '  T<lrl2P*  Ir23l2  *  lr31l2;,  (1.21) 


where  AI^  =  L  -  1^ ,  Equation  (1.23b)  can  be  used  to  write 


AIjAIjAIj 


=  2  r 


12i  ir23l  |r31l  cos(#12  +  *23  +  «3i).  (1.25) 


It  car.  be  seen  from  Equation  (1.25)  that  the  argument  of  the  cosine 
function  involves  the  sums  cc  phase  differences  between  the  three  sam¬ 
pling  points.  Gamo  used  the  argument  of  the  cosine  to  form  a  difference 
equation  through  which  he  expressed  the  phase  of  the  two-point  mutual 
coherence  function.  An  ambiguity  arises,  however,  in  the  sign  of  the 
phase  because  of  the  sign  indeterminacy  of  the  cosine  argument.  In 
effect,  two  images  are  derived,  one  the  reflection  through  the  origin 
of  the  other,  and  supplementary  information  must  be  used  to  resolve  the 
proper  reconstruction. 
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1 . 7  Summary 


We  have  seen  that,  using  the  basic  Huygens-Fresnel  imaging  equation, 
the  fundamental  observable  of  op'.ical  radiation,  the  intensity,  can  be 
used  by  means  of  the  van  Cittert  -  Zernike  equation,  to  infer  the  distri¬ 
bution  of  intensity  over  a  spatially  incoherent  source.  Further,  by 
means  of  the  moment  reduction  formula  of  gaussian  statistics,  we  have 
shown  how  higher-order  moments  can  be  used  to  infer  certain  properties 
of  the  second  moment,  the  mutual  coherence  function.  This  property, 
for  the  fourth-order  moment,  is  inherent  to  the  Hanbury  drown  -  Twiss 
formulation.  Finally,  we  have  seen  the  general  way  in  which  Gamo  used 
the  technique  of  sixth-order  field  correlation  to  derive  the  phase  of 
the  second-order  field  correlation. 

This  introduction  will  serve  as  a  basis  upon  which  we  will  extend 
the  principle  of  intensity  interferometry  to  detection  and  averaging  in 
the  spatial  domain. 


2.  INTENSITY  INTERFEROMETRY  IN  THE  SPATIAL  DOMAIN 
2. 1  Introduction 

In  this  section,  we  will  examine  the  problem  of  using  intensity 
measurements  in  the  far  field  to  infer  the  irradiance  distribution  over 
a  rough  surface.  However,  rather  than  examining  the  two-point ,  time- 
averaged  correlation  of  intensities,  the  far-field  intensity  pattern 
will  be  recorded  for  one  time-resolution  unit  of  the  detector.  The 
spatial  signal  will  then  be  autocorrelated  and  related  to  the 
irradiance. 

We  start  with  a  straightforward  generalization  of  a  method  given 
recently  by  Marchand  and  Wolf.-1  Our  notation  is  similar  and  we  follow 
closely  their  development  through  their  Equation  (35). 


The  Intermediate-Average  Mutual  Coherence  Function 


For  a  stationary  scalar  wave  field,  the  mutual  coherence  function 
for  the  correlation  of  two  space-time  points  is  often  written 


rCij'212*T>T)  5 


2T  j  vT(-r 


t  ♦  T)VT(x2,t)dt, 


(2.1) 


where  x^  =  xni  +  y^j ,  ar.d  the  limits  for  the  time  integration  are 
allowed  to  approach  infinity.  For  this  case,  however,  we  wish  to  keep 
the  parameter  T  finite  and  by  the  subscripts  indicate  that  we  assume  a 
knowledge  of  Vy(x j ,t  +  t)  and  V*(x2,t)  only  over  the  finite  sample  length 

2T.  We  wish  to  call  T(x  ,x.2,t,T)  *he  intermediate-average  mutual 
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coherence  function  and  carefully  stress  that,  for  arbitrary  T  or  shift 
of  origin,  it  may  bear  little  resemblance  to  the  mutual  coherence  func¬ 
tion  defined  by  the  ensemble  average. 

Following  Reference  11,  we  represent  V(x,t)  as  the  temporal  Fourier 
transform  of  the  complex  analytic  signal  (where  a  constant  term  (2rr)~1 
is  suppressed) 


VT(x 


x,t)  =  f  vT(x,u)  exp(-itLt)du),  for  0  <  Jt  |  <  T  (2.2a) 
(3 


-  0  otherwise, 


(2.2b) 


and 


vT(x,w) 


VT (x,t)  exp(ijt)dt. 


(2.2c) 


Substituting  Equation  (2.2a)  into  Equation  (2.1),  interchanging  the 
order  of  integration,  and  time-averaging  and  performing  the  time  inte¬ 
gration,  we  get 


f  (<K.|  »  *2  * 


where 


t,T)  =  ff  wtQL1.£.2>‘Vw2)  exp(-io)1T)sinc[(a)1  -  w2)T]dm1d(i)2, 

°  (2.3) 


_  sin  x 
sine  x  =  — - — 


(2.4) 


and  the  cross- spectral  density  function 


WT(*1>*2><V‘V  5  vTf-l,ul^VT^’“2^ 


(2.5) 


and  the  subscript  T  here  and  later  implies  a  function  based  on  the  elec¬ 
tric  field  statistics  only  for  the  particular  sample  2T  in  length  (here¬ 
after  called  the  detector  resolution  time)  about  the  origin.  The  sine 
function  of  Equation  (2.3)  assumes  the  role  of  a  low-pass  filter.  If  T 
is  very  small,  the  two  frequency  variables  of  Equation  (2.3)  are  essen¬ 
tially  independent  and  all  cross  terms  are  represented  in  the  product  of 
Equation  (2.5).  These  cross  terms  form  a  high-frequency  spectral  con¬ 
tent.  However,  as  T  tends  to  infinity,  the  sine  function  assumes  the 
role  of  a  delta  function,  constraining  correlation  to  occur  only  between 
identical  frequencies  in  the  transform  product  and  forcing  the  integral 
to  a  one-dimensional  form.  In  the  limit  of  large  T,  the  filtered  spec¬ 
trum  of  Equation  (2.5)  becomes  the  mean  square  value  (dc)  of  each  tem¬ 
poral  frequency  component  in  the  signal. 
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Following  Marchand  and  Wolf  and  the  earlier  lead  of  Walther,15 
\T(x,aO  is  represented  in  the  form  of  an  angular  (spatial)  spectrum  of 
plane  waves  in  Cartesian  coordinates  where 


.(x,u)  =  if  aT(p,q,u)  exp[ik(px  +  qy  +  mz)]dpdq  , 

•<x> 

m  =  (1  -  p2  -  q2)*5  if  p2  +  q2  <  1 

=  i  (p 2  ♦  q2  -  l)5*  if  p2  +  q2  >  1, 


k  =  (j/c. 


(2.6) 


(2.7a) 


(2.7b) 


(2.8) 


where  c  is  the  vacuum  velocity  of  light.  Equation  (2.6)  indicates  that 
vT(x,oj)  is  formed  by  a  superposition  of  homogeneous  spatial  waves  propa- 
gating  in  the  half  space  z  >  0  for  the  criterion  expressed  by  Equation 
-and  a  set  evanescent  waves  propagating  parallel  to  the  plane 
z  =  0  for  the  case  described  by  Equation  (2.7b). 

Expressing  Equation  (2.5)  in  the  form  of  the  angular  spectrum  of 
plane  waves  defined  by  Equation  (2.6),  we  get 

03 

*r<i, ‘ffff  VPj.VPz-W"*) 


x  expfikj (PjXj  *  qlXl  +  m^)] 

X  exP[-lk2(P2X2  +  ^2X2  *  m2S2^dPld^idP2dq2’ 


(2.9) 


where 


ATCri-'J,;P2.q2;u,1,m2)  =  a^.q.^a^fp^q  ;a>2).  (2.10) 

Using  Equation  (2.9),  therefore,  the  intermediate-average  mutual 
coherence  function  of  Equation  (2.3)  can  be  written 


r(*l»*2*  t*T)  c  JJ  exp(-iajjT)  sinctCajj  -  wpTjdaijCu. 
0 
to 

"  ffff  Vp.-WW^’ 


C2.ll) 


x  exp[ik1(p1xJ  ♦  q.yj  +  ra^'j] 


x  exp[-ik2(p?x2  +  q2y2  ♦  V2'3dF2dq2dp2dV 


If  T  is  allowed  rc  approach  infinity,  the  limiting  form  of  the  sine  func¬ 
tion  forces  Wj  =  u2j  and  the  intermediate-average  mutual  coherence 
function  clearly  reduces  to  the  form  of  Equation  (2.13)  of  Reference  11 
by  the  elimination  of  one  of  the  time-frequency  integrals. 

2.3  The  Angular  Correlation  Function  and  Other  Spatiai  Correlation 
Functions _ _ _ 

The  cross-spectral  density  function  Wj(x,  ,a)2)  is  now  expressed 

as  a  four-dimensional  spatial  Fourier  integral  as 

QO 

hjCXj , yj ,  Zj  ,x2  ,yz,z2  ;«j  j(u2  )  -  ^ ^ ^  ^T^l  *  gl  '^1  ’  ^2  *  g2  *  ^2  ’*"*1  ,a,2  ^ 


*  txpUCfjXj  +  g^j  ♦  +  g2y2^dfldgldf2ds2*  (2-12) 


Equation  (2.9)  with  Equation  (2.12)  therefore  implies 


V?!  *qi  ;Pz»‘*2J“l»«2)  =  kl2k22,\(klPrkiqi;Zi;“k2P2,"k2q2;Z2;Wr<l32) 


X  expl-HkjiUj  2j  -  k2m2z2)], 

and  specifically  if  Zj  =22=°,  Equation  (2.13)  becomes 


(2.13) 


Ay  (p,  > q j > P2 » q2 1  ~  iO»”^2P2  *"^2q2^'<!,l><1>2^  * 

(2-14) 

Equation  (2.14)  indicates  that  the  angular  correlation  function 


N  I 
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A-j. Cpi»q1;p2*<l2»u,i* ^2^  and  f°ur*c^raensicma*  spatial  Fourier  transform 
of  the  cross -spectral  density  function  are  related  at  the  pl.'ne  z  =  0  if 

fl  =  ki-°i>  f2  =  'k2P2‘  and  g2  =  *M2- 

Further,  the  spatial  transform  of  the  time  spectrum  of  the  field 
can  be  represented  by 

<c 

vT(x,y,z;w)  =JJ  vT(f,g;z;uj)  exp[i(fx  +  gy)]dfdg.  (2.15) 

Comparison  of  Equations  (2.15)  and  (2.6)  indicates 

aT(p,q;w)  =  k' v^kp.kq.z.w)  exp(-ikmz).  (2.16) 


Using  Equation  (2.16)  in  Equation  (2.1)  gives 


WWWV  =  ki2k22  Wi*ki 


x  V-y  ^k2P2  * k2^2  *  22  ’^2  ^ 


(2.1?) 


exp(-i(klmlz1  - 


The  intermediate-average  spatial  correlation  function  is  defined 


VT  -^1  Pi  » ki  3]  »  Z1  '>W1  )V-J  tk2  P2  »k2  q2  *  z2  ’w2  ^ 


=  V^.  (kj  Pj  ,  kj  q^  ,  Zj  ,  k2  p2  ,  k2  q2  ;  ;2  ,u>2  ) 


(2.18a) 


=  ^  P1  *  9j  >  0 ;  k2  p2  ,  k2  q2  ;  0  ;u>j  ,io^ 


expJitkjmjZj  -  k2m2z2)]. 


(2.18b) 


where  the  product  form  of  Equation  (2.18b)  is  implied  by  the  inde¬ 
pendence  of  the  left  side  of  Equation  (2.17)  on  Zj  and  z2.  If  we  set 
Zj  =  z2  =  0,  Equations  (2.17)  and  (2.18)  imply 

2  2 

AjCPj  ,q1;p2,q2;w1  .aij)  =  k1  k2  VT(kJp1  .k^ ;0;k2p2,k2q2;0;u1  ,u>2) .  (2.19) 


Finally,  using  Equations  (2,14)  and  (2.19),  we  find 

vT(fi,g1;0;f2,g2;0;w}»‘*>2)  =  (2. 20) 

Thus  the  relationships  between  the  angular  correlation  function  and  the 
cross-spectral  density  function  are  established  by  Equation  (2.14)  and 
the  angular  correlation  function  and  the  spatial  frequency  correlation 
function  by  Equation  (2.19)  for  the  case  of  intermediate  time  averaging. 

2.4  The  Intermediate-Average  Correlation  Function  in  the  Far  Field 

Now  the  form  of  the  cross-spectral  density. function  is  examined  in 
the  far  field.  Defining  rn  =  (xp2  +  yp2  +  zn2A  we  seek  the  asymptotic 

forms  for  the  case  of  Xj  and  x_0  tending  to  infinity  in  the  paths  indicated 
by  the  direction  cosines 


XJ_ 

r2  ’  r2  ’  r2 


y, 


and 


•1 


r i  r,  rx 

Rewriting  Equation  (2.9)  using  the  definition  of  Equation  (2.10),  we  have 

=  If  aTCpl>qi;wl)  Vl)jdpldqi 

—  an 

QO 

x  JJ  a* (p2,q2;  w2)  expl-ik2(p2x2*  q2y2  +  m222)]dp2dq2.  (2.21) 

As  k,r,  -*■  »  and  fc„r,  ®,  the  asymptotic  form  of  the  two-dimensional 
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integrals  is  given  by  Miyamoto  and  Wolf  as 

2 

!_ 

"2 


V-r~2;  V'V  =  kjk 


f  2tt1 

_  i - L-  C03  a  cos  6. 


/*i  /i  x2  y2  \ 
7T’  r2*  t2’“i’“2 ) 


exptiJkjTj  -  k2r-)J 


(2.22) 


where 


cos  e. ,  — - 

1  r2 


=  cos  e. 


(2.23) 
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and  use  of  Equation  (2.10)  has  been  made.  Finally,  if  Equation  (2.22) 
is  substituted  into  Equation  (2.3),  we  get  for  the  intermediate-average 
mut»  '  coherence  function 


rQil'*2,T,T^  =  fj  e-xP(_iu1T)  sinc((u.  -  a;2)T]da)1dw2 
0 

expfifkjrj  -  k2r2)] 


4it 


ft;  cos  ei cos  e2 


i  2 


rir2 


(2.24) 


A/Xl  X2  >*2  \ 

F7*  fj*  ^  vw2;- 


In  addition,  because  of  the  relationship  given  earlier  relating  the 
angular  correlation  function  to  the  cross-spectral  density  and  the  spatial- 
frequency  functions  [Equations  (2.14)  and  (2.19)],  Equation  (2.221  can 
be  written  in  the  following  forms: 

expfiCk^  -  k 2r 2) ] 

— 


(2.25a) 


2 

(d  )  ~  4tt  cos 

2 

0 j  cos  e2 

Vi  ki>v  „ 

k  x 

2  2 

V 

=  4nz  cos  0j  cos  90  k jk 2 


exP[i (kjr  j  -  k2r2)] 


rlr2 


X 


VT[k 


X1  x2 
__  i*  _•  o-  k  — 

1  r j ’  Ki  tj*  U‘  K2  r2* 


(2.25b) 


Equations  (2.25a)  and  (2.25b)  can  also  be  used  in  Equation  (2.3)  to  pro¬ 
vide  alternate  for>.s  of  Equation  (2.24).  Equations  (2.24)  and  (2.25) 
form  the  modified  version  of  the  forms  given  in  Reference  11,  Equations 
(2.33)  through  (2.35).  Using  these  results,  we  are  now  in  a  position  to 
form  the  self-intensity  function  in  the  far  field. 

2.S  The  Self- Intensity  Function  in  the  Far  Field 


We  now  examine  the  form  of  the  self-intensity  in  the  far  field  by 
letting  points  Xj  =  x_2  =  x^  and  then  letting  the  time  delay,  t,  be 
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zero.  Under  these  conditions,  the  mutual  coherence  function  reduces  to 
the  self-intensity  (Reference  7,  pp.  5Q7-S09)  and,  Rising  Equation  (2.25b) 
in  Equation  (2.3)  and  the  definition  given  In  Equation  (2.18a),  we  have 


I(x,T)  =  4tt^  cos 


eo 

29  J'J'  du)jdu2  — ^ 


i'.w,  exp[i (k.  -  k?)r) 

-!-=•  - i - sinc[  (oij  -  w2)T) 


^j(ki  7*  ki  r:  wl)  ^l(k2  r*  *2  r’  0;  “2)' 


(2.26) 


where,  as  indicated  earlier,  the  sine  function  acts  to  suppress  temporal 
frequencies  in  the  cross  spectrum  higher  than  ^.1/(2T)  Hz.  We  now  utilize 
the  linear  transformation  of  the  time  frequency  variables  (for  which  the 
jacobian  is  unity)  defined  by 


o)j  -  u2  =  0  and 


uij  ♦  u>2 


{2.21) 


Writing  the  u  variables  in  terms  of  these  center-of-mass  coordinates. 


2o  *  e 
2 


2o  -  P 

°2  =  2 


(2.28) 


which,  when  substituted  into  Equation  (2.26),  gives 


2  /*  *°  /*  *° 

I(x,T)  =  cos29  J  da  a2  I  sine  (pT)exp(ipr/c)dp 

O=0  P  =  -*o 

x  i,  Z;0;  0*0/2)  (2-29) 

T\  c  r  c  r  / 

'*/o  -d/2  x  a  -  p/2  y 

-  *T( — —  — —  jbo;  9-o/2j, 


where  the  depx.r.dence  of  the  amplitude  on  the  difference-frequency 
coordinate,  p ,  has  been  dropped,  since  for  quasimonochroraatic  radiation 
o2  >>  j p2/4  J . 

Now,  using  the  defining  transform  relation  of  Equation  (2.15),  we 
write  the  spatial  correlation  function  at  the  source  (z=0)  where 
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0.  co />)  i,  s^l  z-.  o;  a  -  0/2) 


=  — 5-j  /^77"vt(c1,o];0;<!*o/2)  v*(C2,''7;0;o-°/2)  t2-30) 

«  OO 

exp  -i  (  C  (?!  7  +  *i  £)  exp j^i  (^— )  (c2  7  ♦  n2  £)J  d^d^dn. 

Now  the  intermediate-average  spatial-correlation  function, 
v-Uj.n,;  0;  0  ♦  p/2)  v*(Cr.,n^  0;  a-o/2), 


when  considered  with  the  filtering  action  of  the  sine  function  of 
Equation  (2.29)  will  have  an  effective  contribution  only  for  the  low- 
frequency  components  formed  by  the  difference-frequency  terms  M/(2T)  Hz 
or  less. 

in  addition,  we  assume  the  mode  population  to  be  a  slowly  varying 
function  of  a,  since  c  >>  p/2.  This  approximation  can  be  couched  mathe¬ 
matically  by  expressing  the  mode  population  in  a  Taylor-series  expansion 
about  some  center  frequency  ac  and  taking  only  the  first  (constant)  term. 
This  gives  the  idealized  mode  population  the  shape  of  a  rectangle  func¬ 
tion.  Thus  we  write 

V0:„. 0/2)  »;(?,, vO;0-«/2)]1oh  freq 


=  A(t1,n1;o  +  P/2)  A(C2,n2;o-p/2)  exp^^i  -  e2,nj  -  n2;p)]  (2,31a) 


=  AQ^.o)  A(£2,o)  H(p)  expfi^tj  -  i2JP)],  (2,31b) 

where 

II  for  A{£,  c±  p/2)  /  0 
0  otherwise. 

Equation  (2.31)  acknowledges  the  loss  of  the  optical -frequency  phase, 
while  maintaining  the  phase  of  the  intensity  envelope  formed  by  temporal 
beat  modes.  The  degree  to  which  the  phase  of  this  envelope  is  detected 
depends  on  the  bandwidth  of  the  source  and  the  detector  resolution,  2T. 
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Essentially,  these  arguments  were  made  by  Hanbury  3rown  and  Twiss 
(Reference  3b,  p.  311),  except  for  the  defining  of  the  H  function.  Its 
introduction  is  brought  about  by  the  description  of  narrow-band  sources 
by  terms  in  A(£,o).  For  a  thermal  source  of  relatively  large  bandwidth, 
the  maximum  difference  p  will  extend  far  beyond  the  temporal -frequency 
response  of  the  system  (here  reflected  in  the  sine  term  of  Equation  (2.24)] 
and  be  continuous  as  well.  But  for  a  laser  source  exhibiting  a  series  of 
axial  modes,  the  complete  difference-frequency  domain  might  lie  entirely 
within  the  system  response  but  be  piece-wise  continuous  in  its  extent. 

Relative  to  the  representation  of  the  intermediate-average  by  the 
form  of  Equation  (2.31),  we  wish  to  reiterate  a  statement  mude  following 
Equation  (2.1)  that  the  intermediate-averaging  process  may  bear  little 
resemblance  to  the  infinite  time  average,  even  so  far  as  the  detail  of 
the  amplitude  terms,  A (£,o).  This  situation  would  be  serious  if  our 
intent  were  to  infer,  for  example,  the  time- frequency  statistics  of  the 
source,  tx-t  in  the  present  concept,  we  desire  only  to  infer  the  spatial 
properties  of  the  source.  If  we  consider  a  multi -axial -mode  laser  beam 
scattered  from  a  spatially  rough  surface,  the  lack  of  correspondence 
between  the  two  averages  is  unimportant,  for  all  such  mode  history  is 
integrated  out;  all  areas  of  the  scatterer  see  the  sane  mode  characteris¬ 
tics.  Any  mode  fluctuation  would  be  seen  simply  as  a  variation  in 
total  received  power  from  one  sample  to  the  next.  Here,  we  simply 
require  for  one  detector-resolution  time  over  a  spatial  domain  that  the 
process  of  Equation  (2.31)  maintain  the  random  phase  term  $(£j  -  S^.p) 

(due  either  to  the  temporal  mode  structure  of  the  source  or  to  the  scat¬ 
tering  surface  roughness)  with  sufficient  mode  population  [reflected  in 
the  amplitude  terms  A(^, o}]  such  that  quantum  noise  in  both  the  carrier 
wave  and  the  detector  can  be  ignored. 


Using  the  results  of  Equation  (2.31)  in  Equation  (2.29)  and  taking 
6  <<  1,  we  write 


x  exp 


I(P.T)  =  m  c2  do  J*° dp  exp(ipr/c)  sinc(pT)H(p) 

0  -  03 

CD 

■Iff!  A (Ij.o)  AU2,o)  expU^Uj  -£2;d^  (2.32) 

-Co 

[-i  (a^Zi)  («■  ?  *  "l  £)]eX[{i(iT^)(C2r*  n2r)]‘Uld52d,’>d0; 


We  now  have  the  sel^-intensity  in  the  far  field  expressed  as  a  double 
integral  over  sms-  and  difference-frequency  components  as  well  as  two, 
two-dimensional  spatial  Fourier  transforms  over  the  source. 
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Following  Goodman,  '  we  argue  that  the  received  field  at  any  point 
in  the  far  zone  consists  of  a  sum  of  random-amplitude,  random-phase, 
complex  phasors  contributed  by  the  elementary  scatterers.  If  the  size 
of  the  scattering  area  is  large  enough  to  include  many  point  scatterers 
(or  there  are  enough  elementary  coherence  areas  composing  the  source) , 
the  Central  Limit  Theorem  may  be  used  to  conclude  that  the  electric 
field  in  the  detection  plane  is  a  gaussian  random  process  in  a  spatial 
sense. 

Using  the  form  of  Equation  (2.32)  and  its  property  of  spatial  gaus- 
sian  statistics,  we  are  ready  to  form  the  fourth-order  correlation  func¬ 
tion  in  the  far  zone. 

2.6  Fourth-Order  Field  Correlation  in  the  Far  Zone 


In  Section  1.4,  a  theorem  of  fourth-order  gaussian  processes  was 
derived.  This  relation  is  implicit  to  the  work  of  Hanbury  Brown  and 
Twiss,  discussed  in  Section  1.5,  but  was  never  explicitly  utilized.  How¬ 
ever,  KolfS  later  discussed  this  theorem  as  a  plausibility  argument  for 
their  work. 


t  To  form  the  fojrth-order  correlation  function,  we  can  proceed  by 

writing  the  two-point  product  of  intensities  in  the  far  field  using 
Equation  (2.16)  in  a  manner  similar  to  that  of  Hanbury  Brown  and  Twiss. 
However,  to  develop  an  approach  adaptable  to  arbitrary  orders  of  corre¬ 
lations,  as  well  as  to  allow  consideration  of  scattering  surfaces  with 
j  arbitrary  roughness,  we  start  by  writing  the  fourth-order  gaussian 

J  theorem  (here  for  the  electric-field  spatial  variable),  where 

<I(x1)I(x2)>  =  : I (x J ) >  <I(x2)>  ♦  i<V(xi)V*(x2)>!2,  (2.33) 


and  where  the  angle  brackets  <?  indicate  a  spatial  average,  not  the  more 
usual  time  average.  We  can  conclude  from  Equation  (2.33)  that  the 
second-order  intensity  correlation  is  composed  of  two  terms,  of  which 
one  forms  the  square  of  the  second-order  field  correlation;  the  othei 
is  a  spatial  dc  term,  of  no  value  here.  Hanbury  Brown  and  Twiss  elimi¬ 
nated  a  similar  temporal  term  by  means  of  a  dc  block  in  their  electronic 
apparatus.  Part  of  the  temporal  dc  signal  contribution  comes  from  terms 
in  o  =  0.  The  two-point  product  of  intensities  in  the  far-field  (for 
o  =  0)  remains  a  constant  in  the  time  domain  and  is,  therefore,  of  no 
utility  in  the  Hanbury  Brown  -  Twiss  experiment.  This  two-point  product 
does  vary,  however,  as  the  spatial  center-of-mass  coordinate  of  the  two 
points  is  translated  over  the  detection  plane.  Thus,  terms  in  o  =  0 
must  be  evaluated  for  the  case  of  averaging  in  the  spatial  domain. 

Using  the  moment  reduction  properties  of  Equation  (2.33)  along  with 
Equation  (2.32)  and  the  more  general  form  of  the  spatial -frequency  corre¬ 
lation  function  expressed  by  Equation  (2.25b),  we  can  write 


33 


ifr~ TfT-TnMaMftai 
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(2.34a) 


(2.34b) 


exP[-i  (£^)(i1  *£,)]  exp[i  (~^)(l2  •*S8)]dl1dg2|  , 


Next,  we  make  the  following  transformation  to  spatial  center-of- 
mass  coordinates  where 


Is  Li  -  i2  a"d 


(2.35) 


Introducing  these  into  Equation  (2.34),  following  some  algebra,  we  find 
(suppressing  constant  terms), 

Iy*®  i*  ® 

I  a2  da  /  sine (pT)H(o) exp (ipr/c)dp 
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In  going  from  Equation  (2.36a)  to  Equation  (2,36b),  we  have  dropped  the 
Fnurier-transform  terms  in  the  difference  frequency  variable  P ,  since 
they  are  clearly  negligible  in  the  far  field.  The  ensemble  average  has 
been  expressed  as  a  product  of  two  terras*  C(f)  is  a  normalized  phase 
correlation  function'?  describing  the  coherence  interval  over  a  rough 
surface,  and  I(^,w)  is  the  intensity  in  the  global -spatial  variable.  We 
have  made  the  reasonable  assumption  that  the  field  amplitude  is  constant 
-’ithin  a  given  coherence  area  of  the  source;  specifically,  for  £_  suffi¬ 
cient  ly  small  that  C(f_)  ^  0,  A(£  ♦  f/2,o)  =  A(g_  -  f/2,o)  =  A(g^ o)  =  A(£, ct). 
This  final  approximation  is  made  under  the  assumption  that  the  coherence 
area  of  the  source  is  small  relative  to  the  total  source  area,  an  asser¬ 
tion  already  made  in  an  earlier  argument  for  gaussian  statistics.  In 
Equation  (2.36c),  the  circumflex  indicates  a  two-dimensional  spatial- 
Fourier  transform. 


In  order  to  describe  a  spatially  incoherent  surface,  C (f)  is  usually 
allowed  to  assume  the  role  of  a  delta  function. ^  Thus,  Equation  (2.36) 
is  reduced  to  a  single  integral  in  two  space.  Given  this  form,  we  see 
from  Equation  (2.36)  that  the  spatial -averaged,  two-point  intensity  cor¬ 
relation  in  the  far  field  is  proportional  to  the  mod"ius  of  the  spatial 
Fourier  transform  across  a  spatially  rough  surface,  assuming  a  suffi¬ 
ciently  short  exposure  time  2T.  However,  only  for  the  case  that  the 
intensity  distribution  over  the  source  has  even  symmetry  can  the  phase 
of  the  spatial  transform  be  inferred  and  used  to  invert  uniquely  Equation 
(2.36)  to  derive  the  intensity  distribution  on  the  source,  I(£,w). 

Finally,  the  form  of  Equation  (2.36c)  shows  explicitly  that  the 
intensity  autocorrelation  function  in  the  far  field  is  proportional  to 
the  product  of  two  spatial  power  spectra.  The  spatial  power  spectrum 

|i|j  (Xj  -£2)1  j2  is  multiplied  by  |cjj^-  (Xj  +  x2)J  f  ,  the  spatial  power 

spectrum  of  the  correlation  function  describing  the  surface  roughness. 

If  the  surface  is  sufficiently  rough  that  this  function  approximates  a 
delta  function,  tl  en  the  transform  is  essentially  constant,  and  all  spa¬ 
tial  frequencies  o'  the  source  can  be  inferred.  However,  as  tn*.  '-orre- 

lation  interval  increases,  (x^+x^Jf2  acts  to  band  limit  the 

detectable  spatial  spectrum  of  the  source.  This  effect  is  discussed, 
for  example,  by  Kinsly18  for  the  case  of  microdensitometer  imaging  with 
partially  coherent  ligb.‘ 


We  also  wish  to  poir>”  out  that  the  transform  of  the  intensity  dis¬ 
tribution  over  the  sourc* .  lj~  (x^  -x2)J,  is  a  function  of  the  difference 
coordinates  in  the  receive..'  space.  However,  the  transform  of  the  phase 
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correlation  function,  cj^-  (Xj  +  X2)J,  is  a  function  of  the  average 

(global)  coordinates  in  the  receiver  space.  This  terra  cannot  strictly 
come  out  of  the  spatial -averaging  process  of  Equation  (2.36),  since  the 
averaging  is  done  in  the  global  sense.  In  effect,  the  resolving  capacity 
of  a  measurement  in  the  receiver  plane  depends  on  the  position  in  the 
plane.  However,  if  the  coherence  area  in  the  source  plane  is  small,  C 
is  nearly  constant  over  a  large  range  of  its  argument  and,  hence,  is 
generally  insensitive  to  changes  in  the  global  space  variable  at  the 
receiver  around  the  origin.  As  indicated  by  Goodman,17  this  factor  is 
generally  treated  as  a  constant  over  the  receiving  aperture. 

2 An  Illustrative  Experiment 


T?  illustrate  the  above  ideas,  we  suggest  a  simple  experiment 
embodying  these  mathematical  ideas.  A  helium-neon  laser  in  single-axial - 
mode  configuration  is  used  to  transilluminate  a  symmetrical  source 
(  igure  2.1).  The  source  has  a  random-phase  character  to  obtain  spatial 
incoherence.  The  far-tone  intensity  pattern  is  recorded  by  film  using 
an  exposure  time  less  than  the  reciprocal  of  the  source  bandwidth.  This 
exposure  ensures  that  the  spatial  pattern  is  not  washed  out. 


Figure  2.1  A  symmetrical  object  (0)  is  illuminated  with  collimated 
light  from  a  single-axial-mode  laser.  Ground  glass  is  introduced 
at  plane  0  to  achieve  spatial  incoherence.  The  far-field  intensity 
pattern  is  recorded  by  film  at  plane  F.  The  lens  focal  lenqth  is 
indicated  by  f.  3 


Next  the  film  is  developed  so  that  it  is  linear  in  intensity  and  is 
used  to  make  two  identical  positive  transparencies.  The  positives  are 
then  placed  in  a  collimated  beam  (Figure  2.2)  to  form  the  correlated 
intensity  over  an  averaging  area.  The  signal  transmitted  by  the  trans¬ 
parency  pair  is  optically  Fourier-transformed,  a  dc  block  is  inserted 
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Figure  2.2  A  pair  of  transparencies,  Tj  and  T2,  are  made  from  film 
record  in  plane  F  of  Figure  2.1.  Intensity-autocorrelation  function 
is  measured  by  power  meter  for  spatial  lag  f_'  c  Xj  -x^.  DC  block  is 
inserted  between  lenses  at  arrow  to  remove  unwanted  signal. 

to  remove  the  unwanted  average  term,  and  the  total  remaining  irradiance 
is  measured.  This  signal  represents  the  mathematical  expression  given 
by  Equation  (2.36)  for  the  transparency  spatial  lag,  f' '  =  xj  -x^.  Since 
the  source  is  known,  a  priori,  to  be  symmetrical,  the  transform  of  the 
source  intensity  is  pure  real.  The  square  root  of  the  correlation  signal 
(see  Section  3.4  for  details  of  this  operation)  is  proportional  to  the 
spectrum,  which  is  then  known  as  a  function  of  spatial  lag.  Finally, 
this  two-dimensional  signal  is  Fourier-transformed  by  machine  to  give 
the  scaled  source  irradiance. 

2.8  Summary  and  Conclusions 

Having  developed  an  intermediate-average,  mutual -coherence  function 
as  a  starting  point,  we  have  derived  an  expression  for  the  two-point 
intensity  correlation  in  the  far  field,  independent  of  time -aver aging 
except  for  the  temporal  resolution  of  the  detector.  This  result  is  valid 
for  narrow-band,  high-intensity  light  scattered  from  a  spatially  rough 
surface  c:  aibitrary  coherence  area. 

There  are  a  number  of  special  benefits  from  detecting  images  by 
the  technique  of  intensity  correlation.  (1)  The  method^is  relatively 
insensitive  to  the  effects  of  atmospheric  scintillation.-  (2)  Because 
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the  signal  is  detected  in  the  spatial -transform  domain,  high-frequency 
detail  about  the  scattering  surface  translates  to  large  spatial  lags  in 
the  far  field.  This  result  could  be  particularly  important  at  frequen¬ 
cies  where  detector  resolution  is  not  well  developed.  (3)  A  special 
advantage  to  intensity  interferometry  in  the  spatial  domain  is  the  utili¬ 
zation  of  gaussian  statistics  in  the  spatial  (not  temporal)  sense.  By 
this  method,  sources  with  non-gaussian  time  statistics  (such  as  single- 
axial-mode  lasers)  can  be  utilized.  (4)  Still  another  advantage  of 
spatial  detection  is  that  images  of  moving  surfaces  can  be  formed  using 
brief  exposures. 

We  have  therefore  shown  that,  given  a  symmetrical,  spatially  inco¬ 
herent  source  illuminated  by  high-intensity  light,  the  far-zone  intensity 
pattern  can  be  used  to  form  the  optical  image  of  the  source  if  the  signal 
is  recorded  with  sufficiently  short  time  resolution. 


3.  SPATIAL  INTENSITY  INTERFEROMETRY  WITH 
SOURCES  OF  ARBITRARY  SYMMETRY 

3. 1  Introduction 

In  Section  2,  we  showed  that,  given  a  symmetrical,  spatially  inco¬ 
herent  source  illuminated  by  high-intensity  light,  the  far-zone  intensity 
pattern  can  be  used  to  form  the  optical  image  of  the  source  if  the  signal 
is  recorded  with  sufficiently  short  time  resolution.  The  primary  result 
of  this  analysis  is  that  the  far-field  intensity  correlation  function  is 
proportional  to  the  square  of  the  spatial  Fourier  transform  of  the  source 
intensity  distribution.  Sin^e,  in  general,  a  source  exhibits  an  inten¬ 
sity  profile  of  arbitrary  symmetry,  its  spatial  Fourier  transform  is  com¬ 
plex.  The  measurement,  however,  gives  information  only  about  the  spatial 
power  spectrum  (absolute  square  of  the  spatial  Fourier  transform)  and, 
therefore,  only  the  modulus  of  the  Fourier  transform  can  be  inferred. 

The  inversion  of  the  spatial  transform  to  derive  the  source  intensity  is 
thus  impossible  since  the  necessary  phase  information  has  been  lost. 

In  the  experiments  of  Hanbury  Brown  and  Twiss,  the  loss  of  phase  is 
not  a  serious  limitation  since  their  objective  is  simply  the  measurement 
of  star  diameters.  If  a  circular  disk  is  used  as  a  model  for  a  star,  the 
object  is  known,  a  priori }  to  be  symmetrical.  Thus  the  spatial  transform 
of  the  (real)  intensity  is  pure  real.  For  this  situation,  the  phase  of 
the  transform  is  zero  or  n  for  all  spatial  wave  numbers,  and  the  square 
root  of  the  power  spectrum  can  be  taken  (with  a  sign  ambiguity  to  be  dis¬ 
cussed  below)  to  derive  the  spatial  transform  itself. 

In  the  application  of  intensity-interferometric  techniques  to  ter¬ 
restrial  imaging  systems,  though,  the  loss  of  phase  is  a  serious  limi- 

19  20  21 

tation  to  the  method.  A  number  of  authors  *  ’  have  addressed 

themselves  to  the  problem  of  phase  recovery.  Wolf*9  pointed  out  that 
the  complex  transform  of  intensity  has  analytical  properties  that  can 
be  used  in  certain  situations  to  infer  the  phase  from  measurements  of 
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the  modulus  alone.  Gamo  *  has  suggested  a  triple-intensity  correlation 
scheme  by  which  the  phase  can  be  inferred.  However,  there  is  an  ambiguity 
in  the  sign  of  the  phase  angle,  and  two  intensity  profiles  ave  derived, 
one  the  symmetrical  inverse  of  the  other,  and  supplementary  information 
must  be  gathered  to  infer  the  proper-handed  image.  Mehta2-*  has  proposed 
another  scheme  in  which  a  reference  beam  of  an  exactly  known  complex 
degree  of  coherence  is  superimposed  on  the  signal  beam. 

We  suggest  here  a  simple  technique  which,  when  used  with  spatial 
detection,  gives  both  the  amplitude  and  phase  of  the  transform.  This 
method  utilizes  preprocessing  of  the  electric  field  before  detection  to 
exploit  a  symmetry  property  of  Feurier  transform  theory. 

3.2  A  Symmetry  Property 

We  start  by  reviewing  the  well-known  property  that  a  function  ffx) 
can  be  represented  by  a  sum  of  two  functions,24  one  of  which  is  the  even 
(symmetrical)  part  of  f(x),  and  the  other  the  odd.  The  even  [E(x)]  and 
odd  [0(x)]  parts  of  f (x)  can  be  found  simply  by  the  formulas 

E(x)  =  I  [f(x)  +  f(-x)]  (3.1a) 

ard 

0(x)  =  ~  [f(xj  -  f(-x)].  (3.1b) 

If,  in  addition,  we  know  that  ffx)  is  pure  real,  then  the  Fourier  trans¬ 
form  will  be  hermitian;  that  is,  the  transform  of  the  even  part  of  ffx), 

E,  will  be  pure  real,  while  the  transform  of  the  odd  part  will  be  pure 
imaginary.2^ 

By  the  spatial-intensity  method  noted,  the  modulus  of  the  Fourier 
transform  of  the  source  intensity  can  be  derived.  Using  the  above  sym¬ 
metry  properties,  we  preprocess  the  electric  field  before  detection  and 
autocorrelation  so  that  the  source  is  effectively  symmetrized  in  its 
intensity  profile.  Since  its  transform  is  pure  real,  the  power  spectrum 
can  be  used  to  compute  the  spectrum  of  the  even  part  of  the  source  irradi- 
ance.  The  Fourier  spectrum  of  the  even  part  of  the  source  intensity  is, 
however,  the  real  part  of  the  transform  of  the  unpreprocessed  signal. 

Thus,  if  the  modulus  and  the  real  part  of  the  transform  for  each  spatial 
wave  number  are  compared,  the  phase  of  the  transform  can  be  inferred  to 
within  a  sign  and  used  to  compute  the  source  intensity.  The  removal  of 
the  sign  ambiguity  will  be  discussed  below. 

3.3  Mathematical  Model 


We  start  by  writing  a  symmetrized  field  in  the  far  zone  of  the 
source  (£)  plane,  where 

VE(x.)  2  V(x.)  *  V(-Xj).  (3.2) 
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Writing  the  correlation  of  I*(Xj)  and  j'(x2),  where  I '  =  VE(x^)Vg(x^) , 

and  using  the  theorem  of  gaussian  statistics  [Equation  (1.13)],  we  have 


<l'(x1)l'(x2)>  =  i<V£(x1jV*(x2)>r 


(3.3a) 


=  !<V(Xj)V*(x2)  +  V(x1)V*(-x2)  +  Vf-x^V*^)  +  V(-Xj)V*(-x2)>!2  , 

(3.3b) 

where  constant  terms  here  and  in  later  expressions  have  been  suppressed. 
New  the  form  of  Equation  (2.34b)  is  used  to  express  the  intensity  corre¬ 
lation  to  get 

/CO  -  CO 

c2do  J  sinc(pT)H(p)  exp(ipr/c)dp 

tj  —  CO 

00  ( 

*  tiff  |  <A(£rc)A(£2>0)  exp[i<(>(£j  -£2;p)]> 

—  CO 

+  <AUj  ,o)A(-£2,o)  expl-^fCj  ♦  L2’>o)]> 

+  <A(-C1,o)A(C2.o)  exp[i^(-£1  -  C2;p)]>  (3.4) 

♦  <A("£1 » °)f  ('ij.0)  expfi^C-^  +£2;o)]>  j 

«  e*p[-i  -)(«,  'I,)]  «p[i  (£^Zi)  (i2  •i2)j‘*£,«2r. 


We  can  see  now  that  the  form  of  the  intensity  correlation  in  the  far 
field  is  proportional  to  the  square  of  the  spatial  Fourier  transform  of 
fc'jr  terms.  The  nature  of  these  terms  can  be  better  understood  by 
referring  to  Figure  3.1a.  We  show  the  coordinate  axes  for  the  source  (£) 
plane.  Two  points  £,  and  and  their  symmetrical  pair  are  represented. 
The  circle  illustrates  the  coherence  area  characterizing  the  surface 
roughness.  As  discussed  in  Section  2,  the  first  average  of  Equation  (3.4) 
is  zero  unless  points  Cj  -  such  that  they  can  be  enclosed  by  the 
perimeter  of  the  coherence  area.  If  this  criterion  is  satisfied  for  the 
first  average,  it  is  also  met  for  the  last.  Similarly,  the  cross  pro¬ 
ducts  of  Equation  (3.4)  are  zero  unless  =  +£ . 
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(a) 


(b) 


Figure  3.1  Coordinate  axes  in  the  plane  of  the  source,  (a)  Points 
«ind  and  their  symnetrical  pair  are  shown;  center-of-mass  vectors 
are  shown,  where  £  is  the  mean  of  positions  5,  and  ^  and  f  is  the 
difference.  For  the  coherence  area  (CA)  small,  £  =  =  I2  =  L>  and 

f  =  0.  (b)  Source-plane  coordinates  for  =  -j^.  Here  the  center- 

of-mass  vectors  indicate  £  =  0  and  f/2  =  L\  ~  -£.2* 


We  also  wish  to  note  that  the  spatial  average  expressed  by  the 
angle  brackets  on  the  left-hand  side  of  Equation  [3. 4)  can  be  assumed 
to  approach  the  ensemble  average  as  the  area  of  spatial  correlation  ir, 
the  detector  plane  grows  large  with  respect  to  the  correlation  interval 
in  the  same  plane.  The  nature  of  this  correlation  function  is  such  that 
information  about  the  individual  intensity  products  [I' ;T)] 
will  be  averaged  over  the  active  correlation  area.  To  understand  this 
average  more  fully,  we  define  center-of-mass  coordinates  in  the  receiver 
plane  where 


—1  +  —? 

=  - - -  and  f '  =  Xj  -  x_2*  (3.5) 

Thus  the  intensity  correlation  in  the  far  rone  can  be  expressed  func¬ 
tionally  as 

<l'(x1;T)I,(x2;T)>  =  <l'(&‘  ♦  f72;T)l'(g/  -  f’/2;T)>.  (3.6) 

But  for  a  particular  spatial  correlation  lag,  the  difference  between 
receiver  points.  f_  ,  will  remain  constant  while  the  center  of  mass  of 
the  points,  g_' ,  will  vary.  Depending  on  the  domain  over  which  £  is 
varied,  terms  dependent  or.  this  global  variable  should  disappear  leaving 
only  functions  of  the  difference  coordinate. 

We  now  examine  the  first  and  last  terras  (Tl>4)  of  Equation  (3.4). 

If  we  change  the  position  variables  over  the  source,  and  to 
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center-of-mass  coordinates,  £  and  f,  in  the  manner  of  Equation  (3.5), 
we  get  for  these  terms 

CD 

T1>t<  =  ffff  {<A(£  ♦  f/2,o)  A(£  -  f/2,  a)  exp(i$(f)]> 

•  00 

*  <A(-&  -  f/2, c)  A(-g_  ♦  f/2,  o)  exp[i*(-f)]  (3.7) 

x  e>p[-i  ~  (&  ♦  f  /  2 )  •  Xj]  expfi  (g.- f/2)  ‘  x2]dgdf . 


As  discuss*!  in  Section  2,  the  form  of  Equation  (3.7)  implies  that  the 
amplitude  a.iC  phase  of  the  fields  at  the  source  can  be  expressed  in 
product  form.  The  character  of  the  spatial  incoherence  is  determined 
by  the  exponential  terms  within  the  angle  brackets  that  form  the  phase 
correlation  function  over  the  source.  Since  we  assume  a  small  c /..erence 
area  relative  to  the  source  dimensions,  here  in  Equation  (3.7),  £  is 
small  such  that  A(£  +  f/2 ,o)  =  A(£  -  £/2,o)  =  A(j^,o).  Figure  3.1a  illus¬ 
trates  the  center-of-raass  vectors.  For  terms  Tj  4  to  have  value,  £  must 
be  small  so  that  |£j  -  ^ |  is  less  than  the  diameter  of  the  source 
coherence  area.  We  have  also  dropped  the  Fourier  transform  terms  in  the 
temporal  difference  frequency  variable,  p.  With  these  approximations, 
the  first  bracketed  expression  can  be  written 


<A(£  +  £/2»o)  A(£  -  f/2,o)  exp[i<*(£)]> 
=  <(A(£,o)]2>  <exp[i(j> (£)]  > 

=  I(£,«)  C(f), 


(3.8a) 

(3.8b) 


where  the  center-of-mass  temporal  frequency,  a,  has  been  changed  in 
Equation  (3.8b)  to  w.  With  a  similar  development  of  the  second  bracketed 
expression,  Equation  (3.7)  can  be  written 


=  ff  IH&)  *  U-g)3  exp|-i  |  [&  •  (xt 

-  OD 
CO 

x  JJ  C(£)  expj-i  |  f/2  •  (Xj  +xo)|df, 

-OD 

■  ‘e  [r  c[^  (il  * 


(3.9a) 


(3.9b) 
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Here  k  is  the  wave  number  of  the  light,  and  we  have  used  the  evenness  of 
C(f)  to  separate  the  two  integrals.  The  circumflexes  in  Equation  (3.9b) 
indicate  a  spatial  Fourier  transform,  and  IE  *>  I(£)  +  I(-£),  a  symmetrized 
source  intensity  distribution  (since  =  £) . 

In  a  similar  manner,  the  second  and  third  terms  in  Equation  (3.4) 
can  be  examined  to  get 

T2,3  *  iE  [7  <i,  *  i,)]  C  [i  (X,  -  x2)].  (3.10) 

Figure  3.1b  indicates  the  position  of  the  center-of-mass  vectors  for 
this  term.  Here  £j  -  ^  so  that  £  is  small  and  f/2  =  ^  =  £j .  The 

results  of  Equations  (3.9)  and  (3.10)  show  that  the  symmetrization  and 
multiplication  processes  lead  to  two  spatial  Fourier  transforms  of  a 
symmetrized  source  intensity.  One  term  depends  on  the  difference  of 
coordinates  in  the  far-zone  (is  spatially  invariant);  the  other  depends 
on  the  global  term.  Since  the  spatial  average  of  Equations  (3.4)  and 
(3.6)  involves  the  translation  of  the  global  variable  in  the  receiver 
coordinates,  term  T2  3  of  Equation  (3.10)  is  summed  to  a  constant  and 
contributes  only  for’ small  global  sizes  and  can  be  neglected  in  compari¬ 
son  to  T)  4  for  large  global  excursions.  This  result  is  true  only  for 
the  spatial -average  case.  This  can  be  seen  by  noting  that  the  transverse 
receiver  correlation  lag  for  optical  sources  in  the  far  field  is  typi¬ 
cally  a  few  millimetres.  Under  the  assumption  of  a  large  averaging  area, 
this  means  that  term  T2  3  would  contribute  only  for  g'  <  3  mm,  while 
term  T1>4  would  contribute  for  all  up  to  the  limit  imposed  by 

C  [£  (£a  +  **)]■  The  transform  of  the  symmetrized  source  intensity 

expressed  in  Equation  (3.9b)  is  space  invariant.  However,  C,  the  trans¬ 
form  of  the  phase  correlation,  is  a  function  of  the  global  variable.  As 
discussed  earlier,  the  correlation  function  C(f)  approaches  a  delta  func¬ 
tion  if  the  criterion  of  spatial  incoherence  is  applied  to  the  surface 
roughness  over  the  source.  Thus  its  transform,  C(g0,  is  nearly  constant 
over  a  wide  band  of  spatial  frequencies.  Therefore,  C(g()  is  relatively 
insensitive  to  changes  in  its  argument  over  wide  ranges  of  spatial  wave 
number.  But  the  argument  of  C  is,  in  general,  different  for  each  pair 
of  multiplied  intensities  in  the  far  field.  Depending  on  the  bandwidth 
of  C,  the  maximum  averaging  area  can  be  specified  so  that  no  band  limiting 
occurs  for  the  maximum  global  variation. 

Apparently  the  term  of  Equation  (3.10)  arises  because  of  the  way  in 
which  the  source  was  symmetrized.  When  a  source  is  considered  symmetri¬ 
cal,  it  is  the  intensity  only  which  is  symmetrical  about  the  origin. 
However,  here,  by  the  nature  of  the  field  symmetrization  at  the  receiver, 
there  are  two  positions  for  which  the  phase  correlation  function  is  unity 
rather  than  the  usual  one.  The  cross  terms  that  form  T2  3  give  a  scaling 
to  the  Fourier  transform  that  is  quite  different  from  the  term  Tjjl4. 
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Using  the  results  of  Equation  (3,9b)  and  the  fact  that  the  terra  of 
Equation  (3,10)  suras  to  a  negligible  constant  under  a  spatial  average, 
we  can  finally  write  the  correlation  of  intensities  l'(Xj;T)  and  I’fXjjT) 


[I'i.-ipir  i£[£(vs*)]i2’  (3-nb) 


where  Equation  (3.11b)  is  written  with  the  understanding  that 

jC  (Xj  ♦  Xj)J  |2  weights  |iE|2  in  a  non- stationary  way  over  the 

receiver  plane  and  is,  hence,  not  strictly  separable.  However,  this 
term  is  generally  treated  as  a  constant  over  the  receiver  plane.17  In 
this  equation,  IE(£,w)  =  I(£,w)  +  I(-j^,w)  which,  by  Equation  (3.1a),  is 
the  even  part  of  the  source  intensity.  In  Equation  (3.11b), 

|L-  (x  -  x  )  J|2  is  the  spatial  power  spectrum  of  the  even  part  of  the 

2  J  i  -  f  k  q  o 

source  intensity  distribution  and  |C  -^r  (Xj  +  x^)  |  is  the  spatial 

power  spectrum  of  the  source  phase  correlation  function. 

This  information  will  be  used  below  to  form  the  image  I(£,io). 

3.4  The  Transform  Inversion 

Using  the  results  of  Equation  (3.11)  and  the  modulus  of  the  total 
transform,  we  can  now  find  the  intensity  distribution  on  the  source, 
IU,w).  In  Figure  3.2a,  we  illustrate  the  measured  spatial  power  spec¬ 
trum  derived  by  the  method  of  Section  2.  If  the  (positive)  square  root 
is  taken,  the  modulus  of  the  transform  (Figure  3.2b)  is  known. 

Figure  3.3a  illustrates  the  spatial  power  spectrum  of  the  even  part  of 
the  source  intensity  derived  in  Equation  (3.11).  By  the  herraitian 
property  of  the  transform,74  the  spectrum  itself  is  pure  real.  We  take 
the  square  root  (Figure  3.3b),  but  with  a  sign  ambiguity  such  that  two 
spectra  are  derived,  one  the  negative  of  the  other.  The  proper  spectrum 
can  be  inferred,  however,  on  the  basis  of  the  central -ordinate  property 
of  the  Fourier  transform  (Ref.  24,  p.  136),  whereby  the  definite  integral 
of  a  function  over  infinite  limits  in  one  space  is  equal  to  the  value  of 
of  its  transform  at  the  origin.  Since  the  i.,tensity  is  a  positive- 
definite  function,  the  proper  spectrum  is  positive  at  the  origin. 


44 


<  H 


(a)  (b) 

Figure  3.2  (a)  Representation  of  spatial  power  spectrum  of 

source  intensity  distribution,  |i|2,  vs.  spatial  wave  number  kq. 
(b)  Modulus  of  power  spectrum,  jl|,  vs.  spatial  wave  number  kg. 


Figure  3.3  (a)  Pow^r  spectrum  of  even  part  of  source  distribution, 

(b)  Two  possible  real  spectra  of  the  even  part  of  the  source 
distribution.  The  spectrum  which  is  negative  for  k-  =  0  can  be 
eliminated.  b 
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Now,  using  Figures  3,2b  and  3,3b,  we  know  the  modulus  and  the  pro 
jection  on  the  real  axis  for  each  spatial  wave  number  kQg,  Figure  3,4 

Im  I(k0S) 


Re  I(k0S) 


Figure  3.4  For  a  particular  spatial  wave  number,  k05,  the 
modulus  and  its  projection  on  the  real  axis,  !E,  are  shown. 

As  indicated  by  the  lower  dotted  line,  the  modulus  could  be 
located  in  the  lower-half  plane,  thus  giving  two  possible 
values  for  the  phase  of  the  transform. 

illustrates  this  point  by  showing  the  length  of  the  modulus,  [i|,  and  its 
projection  (possibly  negative),  !E,  on  the  real  axis.  It  should  be  noted 
that  there  is  also  a  sign  ambiguity  here  in  the  phase  angle,  for  the 
modulus  could  appear  in  the  upper-  or  lower-half  plane.  From  the  present 
information,  there  appears  to  be  no  way  to  resolve  the  sign  of  the  phase 
angle.  Thus  two  phases  are  derived,  one  the  negative  of  the  other.  If 
each  phase  record  is  used  with  the  modulus  in  a  Fourier  inversion  opera¬ 
tion,  two  images  are  derived,  one  the  symmetrical  inverse  (reflection 
through  the  origin)  of  the  other,  much  the  same  as  Gamo  found  by  a 
triple-correlation  technique. 10> 

However,  in  practice,  the  choice  of  the  proper  (erect)  image  would 
seem  to  be  straightforward  through  the  use  of  a  further  property  of  the 
Fourier  transform.  The  first  moment  of  a  function  can  be  shown  to  be 
(Ref.  24,  p.  138) 

/  xf (x)dx  =  ,  (3.12) 

"'-03 

where  F*(0)  is  the  first  derivative  of  the  transform  of  f(x)  evaluated  at 
the  origin.  It  should  be  noted  that,  her",  origin  refers  to  the  point 
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where  the  argument  of  the  F*  function  is  zero;  it  does  not  imply  that 
tht  derivative  is  inferred  through  a  measurement  constrained  to  the 
neighbornood  of  the  origin  in  the  receiver  coordinates.  Substituting 
I (£)  for  f(x),  the  operation  of  Equation  (3.12)  describes  the  point 
where  1(0  is  mainly  concentrated.  Since,  in  practice,  we  can  always 
translate  the  source  in  a  known  direction  off  axis  by  pointing  the 
receiver  system,  the  source  can  be  moved  into,  say,  the  right-half  plane. 
This  means  that  the  centroid  is  positive,  such  that  the  slope  of  the 
spectrum  must  be  negative  at  the  origin.  Using  this  information  in 
Figure  3,4  for  the  spatial  wave  number  kQg  =  0,  we  can  choose  the  phase 
plot  with  negative  slope.  Thus  the  proper  orientation  of  the  image  is 
inferred,  and  'he  source  intensity  distribution  is  completely  specified. 

We  note  that  in  these  operations  the  Whittaker-Shannon  sampling 
theorem  must  be  considered  with  respect  to  both  the  spatial-frequency 
content  of  the  source  and  its  position  relative  to  the  axis  of  the 
viewing  system.  The  method  here  depends  on  the  uniq  le  inference  of  two 
spectra  by  means  of  the  preprocessed  detection  scheme.  The  behavior  of 
the  spectra  as  inferred  from  the  power  spectra  is  particularly  critical 
as  the  functions  approach  zero  crossings.  There  the  ambiguity  of  the 
functions  must  be  inferred  on  the  basis  of  continuity  arguments.  Thus 
the  limitation  on  the  method  appears  to  depend  on  noise  in  the  process 
which  might  obscure  the  behavior  at  these  critical  points.  Also,  by 
the  effect  of  the  shift  theorem,24  even  a  low  spatial-frequency  source 
when  positioned  off  axis  would  exhibit  fast  phase  variation,  causing 
the  spectra  to  go  through  many  zero  crossings. 

There  is  a  second  method  of  data  reduction  that  is  more  expedient 
than  the  above  dual  operation.  If  the  object  is  pointed  sufficiently 
far  off  axis  such  that  the  twin  images  of  the  symmetrized  field  do  not 
overlap,  then  the  standard  processing  leading  to  the  modulus  of  the 
transform  can  be  ignored.  The  intensity  record  from  the  symmetrized 
field  can  be  used  to  form  the  twin  images,  since  the  transform  is  pure 
real,  using  the  criteria  of  Figure  5.3  to  resolve  the  proper  sign  of  the 
square  root.  After  reconstruction,  one  of  the  images  can  be  discarded. 


3.5  Measurement  Technique 


The  operation  of  Equation  (3.2)  must  be  accomplished  in  real  time 
before  intensity  detection.  Wessely  and  Bolstad,  in  a  study  of 
turbulence- induced  phase  fluctuations,  have  utilized  an  optical  device 
that  achieves  such  field  symmetrization.  A  beam-splitter  cube,  illus¬ 
trated  in  Figure  3.5,  divides  the  incoming  wave  (I)  into  two  parts,  each 
of  which  is  inverted  by  right-angle  prisms  rotated  90*  with  respect  to 
each  other.  After  recombination,  the  field  (at  F)  has  the  form  of 
Equation  (3.2).  We  note  that  the  ray  paths  for  each  leg  of  the  splitter- 
inverter  undergo  an  odd  number  of  reflections;  hence,  with  respect  to 
the  incident  wave,  the  recombined  field  at  F  exhibits  a  change  in  handed¬ 
ness  in  both  the  inverted  and  non-inverted  images.  The  proper  handedness 
can  be  restored,  of  course,  by  utilizing  an  additional  planar  reflection. 
In  an  actual  application,  it  is  likely  that  the  splitter-recombiner  would 
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Figure  3.5  Beam-splitter  cube  that  divides  the 
wave  into  two  parts,  each  of  which  is  inverted 
by  right-angle  prisms  (1,2)  rotated  90°  with 
respect  to  each  other.  After  recombination  at 
F,  the  electric  field  has  the  form  of  Eq.  (3.2). 


be  preceded  by  a  telescope,  a  practical  necessity  to  reduce  unwanted 
background  illumination  and  to  increase  the  effective  size  of  the  detec¬ 
tion  plane.  In  contrast  to  the  experiment  of  Hanbury  Brown  and  Twiss, 
the  optical  parts  preceding  the  detector  plane  would  have  to  be  high 
quality. 

The  inversion  operation  described  in  Section  3.4  can,  of  course,  be 
accomplished  by  electronic  processing  following  reading  of  the  intensity 
records.  It  is  not  clear  how  purely  optical  techniques  could  be  used, 
since  there  is  a  series  of  decisions  that  must  be  made  involving  signs 
of  square  roots  and  continuity  of  functions  at  zero  (abscissa)  crossings. 


Summary  and  Conclusions 


We  have  shown  that  by  preprocessing  the  electric  field  in  a  simple 
way  to  gain  symmetry,  the  technique  of  intensity  interferometry  can  be 
used  to  derive  the  intensity  profile  of  a  source  of  arbitrary  symmetry. 
This  result  is  accomplished  by  comparing  the  intensity-correlation  func¬ 
tions  of  the  processed  and  normal  intensity  records.  It  is  interesting 
to  note  that  it  is  not  necessary  to  invoke  the  methods  of  higher-order 
10  22 

correlations  ’  to  infer  the  phase  of  the  spatial  transform. 
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4.  LASER  SPECKLE  AND  SPATIAL  INTENSITY  INTERFEROMETRY 


4. 1  Introduction 

In  this  section,  will  discuss  the  relationship  between  the  well- 
known  phenomenon  of  laser  speckle  and  spatial  intensity  interferometry. 
There  is  a  strong  simile rity  between  the  results  derived  in  Section  2 
and  previous  analyses  dea’ing  with  laser  speckle;  both  areas  of  theory 
show  that  the  absolute  square  of  the  source-intensity  Fourier  transform 
is  proportional  to  the  far-field  correlation  of  intensities.  Although 
the  mathematical  relations  between  the  source  and  detection  planes  are 
similar,  the  intensity  signals  at  the  detection  plane  that  are  auto- 
correlated  arise  from  different  phenomena.  Laser  speckle  theory  is 
formulated  in  terms  of  a  monochromatic  field  scattered  from  a  rough  sur¬ 
face.  The  time-invariant  speckle  pattern  in  the  far  field  originates 
from  the  interference  of  monochromatic  waves  from  different  points  over 
the  scattering  surface.  Spatial  intensity  interferometry  is  formulated 
in  terms  of  quasimonochromatic  radiation  scattered  from  a  rough  surface. 

A  time-dependent  speckle  pattern  in  the  far  field  is  generated  by  inter¬ 
ference  of  different  temporal  frequencies  scattered  from  the  source. 

Thus  the  temporal  response  of  the  detector  must  be  sufficiently  short  to 
resolve  the  time-varying,  intensity-interferometric  signal.  However, 
the  time- invariant  speckle  pattern  can  be  recorded  with  a  long  exposure. 

Ke  will  discuss  the  mathematical  similarities  between  these  two 
theories. 

4.2  Laser  Speckle  Formulation 

Two  papers,  representative  of  laser-speckle  investigations,  are  due 
to  Goldfischer2^  and  Crane.27  Although  there  are  some  discrepancies  in 
the  results  of  these  efforts,27  the  basic  conclusions  are  similar.  To 
illustrate  the  results  of  laser  speckle,  we  use  the  work  of  Goldfischer. 

We  start  by  noting  that  Goldfischer  argues  for  a  diffuse  surface 
model  that  contains  'an  infinitely  dense  collection  of  scatterers. . . with 
random  phases."  He  h-s  thus  implicitly  argued  for  the  criterion  of  spa¬ 
tial  incoherence  even  though  he  does  not  utilize  this  concept  in  a  quan¬ 
titative  way  [cf.  Equation  (1.7b)].  Goldfischer  begins  with  a  real  form 
of  the  Huygens-Fresnel  equation  (cf.  Section  1.2)  by  writing  the  electric 
field  in  the  receiver  plane  as 

*.»[*(«  *  .  H.D 

where  the  field,  E(x,y),  is  due  to  the  particu  ar  area  on  the  scatterer 
at  a£An,  4^  is  a  random  phase  angle  associated  with  the  scatterer  at 

(£,n),  and  6  is  a  coefficient  written  to  absorb  the  terms  neglected  in 
xy 


Equation  (1.2b)  and  depends  on  the  location  of  the  point  of  observation, 
and  a  is  the  scattering  efficiency  of  the  surface.  We  note  that  Equation 
(4.1)  reveals  a  monochromatic  treatment  of  the  scattered  laser  radiation. 

Next,  the  field  at  the  receiver  plane  is  squared  to  form  the  inten¬ 
sity,  giving 


I(x.y)  =  2^  EL  [l(C»n)I(Ct,n')],J 
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(4.2b) 


The  summations  of  Equation  (4.2)  are  later  allowed  to  approach  integrals 
as  the  areas  A£An  go  to  zero.  There  are  two  distinct  forms  for  this 
equation.  If  the  primed  terms  are  equal  to  the  unprimed  terms,  the  mean 
intensity  in  the  far  field  is  derived.  This  is  the  dc  term  of  no  use  to 
us  here.  The  intensity  I(x,y)  of  utility  here  has  the  mean  subtracted 
out,  and  we  restrict,  as  does  Goldfischer,  the  summations  of  Equation 
(4.2)  to  occur  only  over  the  unequal  primed  and  unprimed  variables,  i.e., 
we  require  C/S'  and/or  n/n'. 

We  now  define  new  variables  of  integration 

~  (C-C')  =  w  and  |  (n-  n’)  =  Jl.  (4.3) 

Using  these  relations,  the  spatial  intensity  autocorrelation  function  in 
the  far  field  of  the  source  is  formed,  where  averaging  eliminates  terms 
involving  random  phases  to  give 


SO 


and  the  spatial  lag  in  the  detector  plane  is  expressed  by 

Y  =  Xj  -  x2  and  6  =  y1  -  y 2.  (4.5) 

Thus  the  right-hand  side  of  Equation  (4.4)  expresses  the  autocorrelation 
of  intensities  in  the  far  field  as  a  Fourier  transform  of  the  auto¬ 
correlation  of  intensities  over  the  source. 

Finally,  using  the  autocorrelation  theorem,  Equation  (4.4)  can  be 
written  in  the  form  given  by  Goldfischer  where 


clQcj,)  I(x,)>  =  (a2/8:rr4)  |i(Y,«)|z,  (4.6) 
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I(Y,*)  =  JJ  d^  dn  I(C.n)  »'xp  [i  ^  (Y£  +  6n)]  •  (4.7) 


This  result,  utilizing  only  a  monochromatic  formulation,  is  identical  in 
form  to  the  comparative  expression  of  spatial  intensity  interferometry 
given  in  Equation  (2.36);  the  intensity  autocorrelation  function  in  the 
far  zone  of  a  spatially  incoherent  source  is  proportional  to  the  absolute 
square  of  the  spatial  Fourier  transform  over  the  source.  However, 

Equation  (4.6)  reflects  no  constraint  involving  the  exposure  time.  As 
indicated  by  Equation  (4.2b),  the  useful  intensity  signals  in  the  receiver 
plane  arise  from  interference  of  monochromatic  waves  over  different  por¬ 
tions  of  the  source  plane. 


Summary  and  Conclusions 


We  have  seen  that  both  laser  speckle  and  spatial  intensity  inter¬ 
ferometry  techniques  indicate  that  the  far-field  correlation  of  intensi¬ 
ties  is  proportional  to  the  absolute  square  of  the  source-intensity 
Fourier  transform.  However,  laser  speckle  theory  describes  scattered 
monochromatic  radiation  in  a  time- independent  pattern  requiring  no  expo¬ 
sure  control  to  detect.  In  contrast,  spatial  intensity  interferometry 
is  formulated  in  terms  of  quasimonochromatic  radiation  that  yields  a 
time-varying  pattern.  Exposure  control  commensurate  with  the  bandwidth 
of  the  radiation  is  required  to  resolve  the  pattern. 


We  have  seen  that  both  the  time-varying  and  time- invariant  speckle 
patterns  contribute  to  the  useful  signal  describing  the  source  intensity 
distribution.  We  must  ask,  is  it  to  be  expected  that  speckle  recording 
with  short  exposures  would  result  in  an  increased  signal  in  the  auto¬ 
correlation  process  due  to  the  time-varying  signal  which  would  otherwise 
be  lost?  It  would  seem  that  the  answer  lies  in  the  nature  of  the  corre¬ 
lation  between  the  two  kinds  of  patterns.  This  question  will  be  examined 
experimentally  in  the  next  section. 


5.  EXPERIMENTS  IN  INTENSITY  CORRELATIONS 
S.l  Introduction 

In  this  section,  we  examine  supporting  experimental  evidence  for 
some  of  the  basic  ideas  encountered  in  this  dissertation.  Specifically, 
in  Section  2,  we  showed  that  through  the  utilization  of  intensity  beat 
signals  at  a  detection  plane,  a  spatial  Fourier  transform  relation  is 
derived  with  respect  to  the  source-plane  irradiance  distribution.  In 
Section  4,  we  showed  that  the  Fourier  transform  relation  between  the 
source  and  detection  planes  using  time-varying  intensity  beats  of  an 
optical  source  is  identical  to  that  derived  in  laser  speckle  theory2***2? 
using  a  monochromatic  formulation.  Since  these  two  theories  are  identi¬ 
cal  in  their  spatial  transform  relations,  we  will  show  first  the  results 
of  unpublished  laser  speckle  experiments  performed  by  Peppers28  which 
verify  these  relations. 

As  we  have  discussed  in  Section  4,  the  signal  utilized  in  laser 
speckle  theory  (and  by  Peppers  in  his  experiments)  is  time  invariant, 
and,  hence,  can  be  recorded  with  a  long  exposure.  Here,  we  will  refer 
to  this  signal  as  the  "dc  speckle  pattern."  By  analogy,  the  time-varying 
intensity  signal  of  import  to  the  intensity- interferometry  formalism  we 
will  refer  to  as  the  "ac  speckle  pattern."  In  the  latter  half  of  this 
section,  we  will  describe  our  own  measurements  of  intensity  fluctuations 
in  a  laser  speckle  pattern.  By  these  experiments,  we  can  infer  the 
existence  of  an  ac  speckle  pattern,  the  nature  of  its  fluctuations,  and 
its  relation  to  the  dc  speckle  pattern. 


5.2  Experiments  in  Spatial  Intensity  Correlation 


28 

We  describe  now  the  intensity  correlation  experiments  cf  Peppers 
which  illustrate  the  transform  relations  given  in  Sections  2  and  4.  The 
experimental  configuration  that  he  used  is  similar  to  that  illustrated 
earlier  in  Figure  2.1.  A  helium-neon  laser  has  been  used  to  illuminate 
an  opal-glass  diffuser  (to  introduce  spatial  incoherence)  followed 
immediately  by  an  object;  it  consists  of  a  200-mesh  grid  of  square  holes 
outlined  by  a  1.22-nm  circular  hole  in  a  metal  plate.  A  photograph  of 
this  object  is  shown  in  Figure  5.1.  This  grid,  along  with  the  glass 
diffuser,  was  placed  at  plane  0  illustrated  in  Figure  2.1.  At  plane  F, 
the  resulting  dc  speckle  pattern  was  recorded  photographically  and  is 
shown  here  in  Figure  5.2. 
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The  rnotographic  record  made  at  plane  F,  Figure  2.1,  was  used  to 
make  a  pa;r  of  positive  transparencies  in  order  to  form  the  autocorre¬ 
lation  function  given  in  Equation  (4.6).  In  Figure  2.2,  we  illustrated 
a  standard  way  in  which  an  autocorrelation  function  can  be  formed  by 
translating  physically  one  transparency  to  all  vector  spatial  lags  f ' 1 . 
An  alternative  approach,  used  by  Peppers,  is  shown  in  Figure  5.3. 
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Figure  5.3  Optical  system  used  to  form  the  correlation  of 
transparencies,  Ti  and  T2.  Collimated  light  impinges  from 
the  left  onto  diffuser,  D.  Distance  between  Tj  and  T2  is  d, 
and  f  indicates  the  focal  length  of  the  lens.  *F  denotes 
the  film  plane. 

Collimated  light  from  the  left  is  diffused  at  plane  D  where  the  first 
transparency,  Tj,  is  placed.  At  a  distance  d  from  the  first  transparency, 
a  second  identical  positive,  T2,  is  positioned  and  followed  directly  by 
a  lens  of  focal  length  f.  At  the  focus  of  the  lens,  film  (FI  records 
the  two-dimensional  autocorrelation  function.  Reflection  on  the  geometri¬ 
cal  optics  interpretation  of  this  system  reveals  this  operation  as  well 
as  the  equation  governing  the  spatial  magnification;  in  the  film  plane, 

s=ff’/d,  (5..) 

where  £  is  the  spatial  lag  in  the  film  plane,  f  is  the  focal  length  of 
the  lens,  f’  is  the  spatial  lag  between  the  transparencies,  Tj  and  T2, 
and  d  is  the  distance  between  the  two  transparencies. 

The  results  of  this  operation  are  shown  in  Figure  5.4  We  note  a 
central  spot  containing  low-frequency  information  about  the  object  cen¬ 
tered  among  four  side  lobes,  reflecting  sv.rong  spatial-frequency  content 
due  to  the  periodicity  of  the  200-mesh  grid.  The  signal  illustrated  in 
Figure  5.4  is  the  function  given  by  Equation  (4.6).  It  should  be  noted 
that  the  intensity  relation  of  Equation  (4.6)  is  similar  to  the  expres¬ 
sion  in  electric  fields  given  by  Equation  (1.3).  Peppers  recorded  the 
square  of  the  electric  field  expressed  by  Equation  (1.3)  by  using  the 
Courier  transform  apparatus  of  the  kind  illustrated  in  Figure  2.1  after 
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Figure  5.4  Optical  autocorrelation  of  speckle  pattern  shown  in  Figure  5.2 

first  removing  the  ground  glass  at  plane  0.  The  result  is  presented  in 
Figure  5.5  and  is  to  be  compared  with  the  function  shjwn  m  Figure  5.4. 
Except  for  the  spatial  noise  of  Figure  5.4,  the  two  functions  are  similar 
in  the  manner  indicated  by  Equations  (1-3)  and  (4.6). 


Figure  5.5  Diffraction  pattern  of  object  in  Figure  5.1  made  by 
configuration  illustrated  in  Figure  2.1  after 
removal  of  optical  diffuser  at  plane  0.^8 


28 
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In  the  previous  chapter,  the  autocorrelation  theorem  was  used  by 
Goldfischer26  to  relate  the  Fourier  transform  of  the  source  intensity 
autocorrelation  [Equation  (4.4)]  to  the  equivalent  expression  of  the 
squared  modulus  of  the  source  intensity  Fourier  transform  [Equation  (4.6)]. 
Peppers  illustrated  this  equivalence  by  the  following  operations.  First, 
the  source  irradiance  (Figure  5.1)  was  photographea  nd  autocorrelated 
to  form  the  object  autocorrelation  function  (Figure  S.b).  Next,  both 
the  intensity  autocorrelation  function  (Figure  5.4)  and  the  squared  modu¬ 
lus  of  the  diffraction  pattern  (Figure  5.5)  were  optically  transformed 
to  effect  the  Fourier  inverse  of  the  first  operation  of  Equation  (4.4); 
these  results  are  shown  in  Figures  5.7  and  5.8,  respectively.  Although 
the  signal  quality  of  Figure  5.7  is  somewhat  poor,  the  similarity  among 
these  operations  can  be  seen. 

On  the  basis  of  these  intensity  correlation  experiments,  it  appears 
that  the  spatial  transform  relations  predicted  by  Goldfischer  for  the  dc 
speckle  pattern  can  be  verified.  We  have  shown  in  this  dissertation  the 
equivalence  between  the  spatial  transform  relations  for  the  dc  and  ac 
speckle  signals.  It  remains  now  to  show  the  specific  relationship 
between  these  two  speckle  patterns.  We  examine  this  subject  experimen¬ 
tally  in  the  next  sections. 

5.3  Conditions  for  Coherence 

Before  examining  laser  speckle  patterns  experimentally,  we  find  it 
useful  to  review  the  conditions  for  obtaining  spatially  fixed  inter¬ 
ference  patterns.  A  discussion  of  these  criteria,  as  well  as  detailed 
experiments  into  the  nature  of  speckle  patterns,  has  been  given  by 
Martienssen  and  Spiller.^9  We  follow  closely  their  development. 

OQ 

It  can  be  asserted  that  each  oscillation  mode  of  a  source  can 
interfere  only  with  itself.  Since  the  volume  of  an  oscillation  mode  in 
phase  space  is  h3  (where  h  is  Planck’s  constant),  the  Heisenberg  uncer¬ 
tainty  relation  can  be  used  to  write  the  ranges  in  position  (x,y,z)  and 
momentum  (px,p.,p?)  for  a  photon  in  the  mode,  where 


Ax  Apx  =  h. 

(5.2a) 

A.V  Apy  =  h. 

(5.2b) 

Az  Apz  =  h. 

(5.2c) 

If  the  mode  propagates  in  the  z  direction,  the  uncertainty  of  the  mcmen- 
tura  in  the  x-z  and  y-z  planes  [per  Equations  (5.2a)  and  (5.2b)]  can  be 
expressed  in  terms  of  the  angular  apertures  2Aux,  2Auy  of  the  beam  in 
the  x-z  and  y-z  planes.  Using  the  relations  h/p  =  X  «nd  sinAux  =  Apx/p, 
sinAUy  *  Apy/p,  we  can  write 


Figure  5.6  Autocorrelation  of  object 
function  illustrated  in  Figure  5. I.28 
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Figure  5.7  Optical  Fourier  transform 
of  signal  function  of  Figure  5. 4. 28 


Figure  5.8  Optical  Fourier  transform 
of  signal  function  of  Figure  5. 5. 28 


Ax  sin  Au  =  A, 
x 

(5.3a) 

Ay  sin  Au  =  A. 
y 

(5.3b) 

Using  Ap^,  =  hAv/c  and  Az  =  cAt,  where  Av  is  the 
light  and  At  is  the  time  during  which  the  beam 
(5.2c)  can  be  used  to  get 

frequency  spread  of  the 
is  observed.  Equation 

Av  At  =  1 . 

(5.4) 

The  values  of  the  angular  apertures  2Au  ,  2Au  which  fulfill  Equations 
(5.3a)  and  (5.3b)  are  called  the  coherence  angles,  and  the  value  of 
observation  time  At  fulfilling  the  criterion  of  Equation  (5.4)  is  called 
the  coherence  time.  Angles  less  than  the  angular  apertures  define  the 
conditions  for  spatial  coherence,  or 

a  <<  Au  , 
x  x 

(5.5a) 

and 

a  <<  Au  . 

y  y 

(5.5b) 

Observation  times  less  than  the  coherence  time 
for  temporal  coherence,  or 

At  fulfill  the  condition 

T  <<  At. 

(5.6) 

If  the  criteria  of  Equations  (,5.5a)  and  (5.5b)  are  met,  a  section  of  a 
light  beam  will  experience  no  fluctuation  in  space,  although,  depending 
on  the  observation  time  T  and  the  inequality  expressed  by  Equation  (5.6), 
it  may  experience  fluctuations  in  time.  In  addition,  if  a  light  beam  is 
observed  in  a  time  fulfilling  the  criterion  of  Equation  (5.6)  (is  tem¬ 
porally  coherent),  then  there  will  be  no  fluctuations  in  time;  but  if 
the  observation  area  is  large  enough,  there  will  be  fluctuations  in  space 

5.4  Experiments  in  Temporal  Intensity  Fluctuations 

29 

Martienssen  and  Spiller  were  apparently  the  first  to  suggest  the 
use  of  a  laser-illuminated  rotating  ground  glass  to  model  a  thermal  source 
with  variable  coherence  time.  Even  with  a  laser,  the  finite  width  of  its 
frequency  spectrum  implies  some  fluctuations  in  intensity.  These  fluc¬ 
tuations  can  be  filtered  using  a  detector  with  a  slow  response.  If  the 
laser  illuminates  a  section  of  ground  glass,  the  random  phase  variations 
introduced  in  the  beam  are  similar  to  the  phase  mapping  that  would  be 
observed  ever  an  incoherent  source  if  the  measurement  were  made  in  a 
time  less  than  the  coherence  time.  Thus  the  stationary  ground  glass 
corresponds  to  the  case  of  infinite  temporal  coherence.  By  rotating  the 
ground  glass,  the  phase  variations  over  the  source  can  be  made  to  change 
with  a  coherence  time  inversely  proportional  to  the  angular  velocity  of 
the  ground  glass. 
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This  experimental  approach  was  used  by  Martienssen  and  Spiller  to 
repeat  the  Hanbury  Brown  -  Twiss  laboratory  experiment  in  which  the  cor¬ 
relation  of  intensity  fluctuations  in  a  beam  was  examined  as  a  function 
of  detector  separation  [see  Equation  (1.21)].  We  have  duplicated  that 
experiment  using  a  setup  illustrated  in  Figure  5.9.  A  Spectra-Physies 
Model  125  helium-neon  laser  is  used  to  illuminate  a  section  of  ground 
glass  that  can  be  kept  stationary  or  rotated  by  a  motor.  The  far-field 
pattern  is  divided  by  means  of  a  beam  splitter  and  examined  by  two 
photomultipliers,  one  of  which  can  be  translated  transverse  to  the  opti¬ 
cal  axis.  The  outputs  of  the  photomultipliers  are  fed  through  ampli¬ 
fiers  with  a  dc  block  into  an  electionic  correlator  (Princeton  Applied 
Research  Correlator,  Model  101).  The  cross-correlation  of  the  two  sig¬ 
nals  is  observed  by  means  of  an  oscilloscope  readout. 


Figure  5.9  Experimental  configuration  for  the  correlation  of  temporal 
fluctuations  within  a  laser  speckle  pattern.  Ground  glass  (GG)  can  be 
spun  by  means  of  motor  (M) ;  beam  splitter  (S)  divides  speckle  pattern 
into  two  images  that  are  viewed  by  photomultipliers  (PM1.PM2);  outputs 
are  fed  through  amplifiers  (A1,A2)  into  electronic  correlator  and  into 
display  oscilloscope.  Photomultiplier  PM1  can  be  translated  across 
the  beam. 

In  Figure  5.10,  we  show  the  correlator  output  for  the  case  of  the 
rotating  ground  glass  and  complete  alignment  of  the  two  photomultipliers. 
The  correlation  of  electronic  signals  is  shown  as  a  function  of  time 
delay  between  the  two  signals.  The  integration  time  of  the  correlator 
is  iS  s,  and  the  maximum  delay  shown  is  50  ms.  The  computed  correlation 
time  of  the  irradiance  fluctuations  can  be  shown  to  be  on  the  order  of 
20  us.  Thus  the  correlation  behavior  of  central  interest  is  not  resolved 
by  the  scale  used  in  Figure  5.10.  Because  of  the  rotation  of  the  ground 
glass,  the  correlation  function  is  periodic.  This  periodicity,  as  well 
as  the  fine  structure  on  the  correlation  function  due  to  the  non-ideal 
operation  of  the  spun  ground  glass,  can  oc  easily  removed  by  use  of 
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Figure  5-10  Oscilloscope  trace  of  correlator  output.  Corre¬ 
lation  C(t)  is  plotted  vs.  time  delay,  i,  between  fluctuating 
signal  components.  Both  optical  detectors  viewed  the  same 
area  of  the  laser  speckle  pattern.  Each  division  along  the 
t  axis  corresponds  to  a  time  delay  of  10  ms. 

high-pass  filters  before  the  process  of  correlation.  In  the  experiments 
that  follow,  the  correlation  lag  (t)  scale  was  expanded  to  give  the 
maximum  (system  tested)  temporal  resolution  on  the  order  of  5  us.  In 
addition,  the  following  measurements  make  use  of  the  correlation  at  zero 
time  delay  (origin  of  the  correlation  function). 

We  examined  the  correlation  of  intensity  fluctuations  as  a  function 
of  detector  separation.  Figure  5.11  shows  the  results;  the  correlation 
is  given  versus  spatial  lag  (Jf').  The  spatial  correlation  interval  in 
the  receiver  plane  is  a  function  of  the  spot  size  on  the  ground  glass; 
for  the  setup  used  the  predicted  interval  is  about  0.065  mm.  The  detector 
apertures  were  about  0.3  mm  for  reasons  of  detector  efficiency  in  ensuing 
experiments.  Thus  the  shape  of  che  correlation  is  more  a  function  of 
detector  aperture  than  source  function.  The  apparent  correlation  interval 
is  about  0.3  mm.  It  can  be  seen  that  the  results  show  a  residual  corre¬ 
lation  level  which  is  due  to  low-frequency  components  in  the  beam.  These 
low  frequencies  were  progressively  eliminated  using  the  high-pass  capa¬ 
bility  of  the  amplifiers.  The  residual  correlation  was  essentially  gone 
after  removal  of  fluctuations  lower  than  500  Hz.  These  results  illustrate 
the  predicted  relation  of  Equation  (1.21)  and  compare  as  well  with  the 
experiment  of  Martienssen  and  Spiller  (Ref.  29,  Figure  7). 

The  question  posed  in  an  earlier  section  involves  the  lelation 
between  the  static  (dc)  and  time- varying  (ac)  speckle  patterns  due  to  a 
laser  source.  The  setup  of  Figure  5.9  was  used  to  investigate  this  phe¬ 
nomenon  by  keeping  the  ground  glass  stationary  and  examining  the  temporal 
fluctuations  due  to  the  laser  itself.  The  stability  of  this  helium-neon 
laser  is  such  that  the  intensity  fluctuations  compose  only  a  small 
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Figure  5.11  Normalized  crosscorrelation  of  intensity  fluctuations 
vs.  difference  of  detector  positions  (_f * )  for  the  case  of  spinning 
ground  glass.  Spatial  correlation  interval  (Cl)  and  detector 
aperture  size  (AS)  are  shown.  Coherence  time  of  the  radiation 
^20  ys.  Curves  a,  b,  and  c  indicate  high-pass  filtering  starting 
at  1,  100,  and  500  Hz,  respectively. 

fraction  (<  1%)  of  the  average  intensity  lev?l.  One  might  imagine  that 
the  temporal  mode  fluctuations  across  the  laser  wave  front  are  every¬ 
where  identical  since  the  source  is  spatially  coherent.  In  the  far 
field  where  tne  speckle  pattern  appears,  the  temporal  intensity  fluctu¬ 
ations  should  be  completely  correlated;  this  concept  can  be  stated 
simply 

I(Xj,t)  =  A(x2)  I(x2,t),  (5.7) 

or  the  temporal  fluctuations  of  intensity  at  one  point  in  the  speckle 
pattern  are  identical  to  the  temporal  fluctuations  at  another  point  in 
the  pattern  times  a  constant  which  depends  on  the  average  value  of  the 
intensity  at  that  point.  This  constant  is,  in  fact,  the  dc  speckle 
pattern  at  that  particular  point. 

The  crosscorrelation  function  was  measured  as  before  but  with  the 
ground  glass  stationary.  For  each  spatial  lag,  the  correlation  at  zero 
temporal  lag  and  the  dc  photocurrent  on  the  movable  detector  were  meas¬ 
ured.  We  note  that  the  intensity  fluctuations  due  to  the  laser  are  a 
full  two  orders  of  magnitude  lower  than  the  rotating  ground  glass  fluc¬ 
tuations.  The  photomultipliers  were  operating  near  the  limit  imposed 
by  shot  noise;  evidence  for  this  observation  is  given  by  comparison  of 
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the  crosscorrelation  of  photomultiplier  signals  at  zero  spatial  lag  with 
the  autocorrelation  of  one  signal.  The  autocorrelation  is  a  factor  of 
two  greater  because  of  the  correlation  of  noise;  under  crosscorrelation 
the  uncorrelated  noise  terms  vanish.  The  correlation  was  normalized  by 
the  dc  value,  and  the  results  are  shown  in  Figure  5.12.  As  before, 


Figure  5.12  Normalized  crosscorrelation  of  intensity  fluc¬ 
tuations  vs.  difference  of  detector  positions  (Jf ' )  for  the 
case  of  static  ground  glass.  Spatial  correlation  interval 
and  detector  aperture  size  are  shown  in  Figure  5.11. 

Curves  a  and  b  indicate  high-pass  filtering  of  500  and  1  Hz, 
respectively. 

high-pass  filtering  of  1  and  500  Hz  was  applied  to  the  two  data  runs 
shown.  The  detector  aperture  size  and  spatial  correlation  interval 
remained  the  same.  As  can  be  seen,  the  temporal  fluctuations  are  com¬ 
pletely  correlated  for  distances  greater  than  forty  times  the  spatial 
correlation  interval  for  this  situation  in  which  a  spatially  coherent 
source  illuminr.tes  a  time-independent  phase  screen.  Here,  the  ac  speckle 
pattern  is  identical  to  the  dc  speckle  pattern  and  appears  merely  as  a 
fluctuation  around  the  average  (dc)  value.  On  the  average,  then,  in 
this  situation  there  is  no  gain  by  making  a  fast  exposure,  even  if  the 
depth  of  modulation  were  somewhat  greater  than  one  percent. 

It  should  be  noted  that  this  condition  was  anticipated  earlier  above 
Equation  (2.32).  There  it  was  argued  that  all  areas  of  the  scatterer 
see  the  same  mode  history  (i.e.,  the  source  is  spatially  coherent)  and 
that,  therefore,  the  mode  fluctuation  of  the  laser  would  simply  be  seen 
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as  a  variation  in  total  received  power  from  one  sample  to  the  next.  The 
measurements  of  Figure  5.12  bear  out  this  concept. 

5.5  Summary  and  Conclusions 


In  this  chapter,  we  have  seen  the  way  in  which  the  far-field  speckle 
pattern  of  a  spa*ially  incoherent  source  can  be  related  to  the  intensity 
over  the  source  itself.  We  have  also  discussed  the  conditions  for  spatial 
and  temporal  coherence  and  investigated,  experimentally,  the  correlation 
of  temporal  fluctuations^in  laser  speckle  patterns.  It  can  be  seen  that 
laser  speckle  theory26,2/  forms  the  limiting  case  in  intensity  inter¬ 
ferometry  in  which  the  coherence  time  of  the  source  becomes  large  and 
for  which  the  averaging  is  accomplished  in  a  spatial  sense.  At  the  oppo¬ 
site  extreme  is  found  the  Hanbury  Brown  -  Twiss  experiment  in  which  an 
incoherent  source  is  monitored  by  means  of  temporal  fluctuations  averaged 
along  the  time  axis.  We  have  seen  [in  Equation  (1.20)]  that  many  of  the 
beat  frequencies  in  the  carrier  wave  are  not  useful,  practically,  because 
of  the  limitations  in  detector  frequency  response.  This  dissertation 
bridges  these  limiting  cases  by  formulating  the  case  of  quasimonochromatic 
radiation  (finite  coherence  time)  detected  and  averaged  in  the  spatial 
domain. 

To  elucidate  the  ultimate  equivalence  among  these  approaches,  we 
turn  to  Figure  5.13.  We  have  represented  an  ensemble  of  similar  con¬ 
figurations  of  source-detector  planes.  We  imagine  the  field-amplitude 
distribution  across  each  plane  to  be  identical  as  well  as  the  two  detec¬ 
tion  points  in  the  receiver  planes.  The  Hanbury  Brown  -  Twiss  case 
(quasimonochromatic  source,  temporal  averaging)  is  taken  first.  We  let 
each  member  of  the  ensemble  represent  an  intensity  measurement  made  in  a 
time  less  than  the  coherence  time  of  the  radiation.  The  phase  structure 
over  each  source  surface  is  therefore  frozen,  but  statistically  inde¬ 
pendent  from  one  another.  Hanbury  Brown  and  Twiss  started  by  forming 
the  product  of  intensities  at  the  two  detection  points  in,  say,  the  ith 
member.  They,  of  course,  wanted  the  ensemble  average  of  that  product 
and,  through  the  assumption  of  temporal  stationarity  and  ergodicity, 
sought  the  equivalent  operation  by  means  of  a  time  average.  This  can  be 
visualized  by  allowing  the  member  j  to  represent  the  identical  source 
amplitude  now  exhibiting  a  new  phase  structure  due  to  the  evolving  tem¬ 
poral  mode  histroy  over  the  source.  As  many  samples  are  averaged  in 
time,  the  function  approaches  the  ensemble  average. 

In  the  case  of  spatial  averaging,  the  ith  sample  can  be  taken  as  a 
spatial  recording  of  intensities  in  the  x-plane  over  a  time  short  com¬ 
pared  with  the  coherence  time  of  the  radiation.  Two  points  within  the 
pattern  can  be  sampled  to  form  the  product  but,  again,  it  is  the  ensemble 
average  that  is  desired.  Here,  through  the  assumption  of  spatial  sta¬ 
tionarity  (and  a  sufficiently  large  area  of  spatial  averaging),  the  spa¬ 
tial  average  is  assumed  equivalent  to  the  ensemble  average. 


Figure  5.13  An  ensemble  of  similar  experiments.  Source  ampli¬ 
tude  and  receiver  points  are  identical.  Random  phase  variations 
over  each  source  plane  are  statistically  independent. 


We  note  that  all  the  schemes  involving  intensity  correlation  require 
a  random  phase  structure  across  the  source.  It  is  achieved  naturally 
for  an  incoherent  source,  and  for  a  spatially  coherent  laser  it  can  be 
accomplished  by  use  of  a  diffusing  glass.  We  note  that  this  condition 
is  described  analytically  by  the  exponential  terra  of  Equation  (2.31). 

This  term  is  written  sufficiently  general  to  include  phase  perturbation 
due  either  to  spatial  surface  roughness  or  to  temporal  mode  difference. 

It  should  be  noted  that  a  suitably  prepared  broad-band  laser  might  be 
utilized  without  ground  glass  if  the  temporal  mode  structure  could  be 
sufficiently  varied  in  a  statistically  stationary  way  over  the  scattering 
surface. 

We  note  again,  as  we  did  at  the  beginning  of  Section  1,  that  the 
focus  of  primary  interest  in  optical  imaging  is  the  spatial  domain.  How¬ 
ever,  we  have  seen  that  the  ability  to  detect  intensity  fluctuations  and 
the  manner  in  which  an  average  is  finally  taken  rests  in  a  most  crucial 
way  with  the  nature  of  the  temporal  statistics  of  the  radiation. 
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